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Abstract 


In the spectral stability analysis of localized patterns to singular perturbed evolution prob¬ 
lems, one often encounters that the Evans function respects the scale separation. In such cases 
the Evans function of the full linear stability problem can be approximated by a product of 
a slow and a fast reduced Evans function, which correspond to properly scaled slow and fast 
singular limit problems. This feature has been used in several spectral stability analyses in or¬ 
der to reduce the complexity of the linear stability problem. In these studies the factorization 
of the Evans function was established via geometric arguments that need to be customized 
for the specific equations and solutions under consideration. In this paper we develop an al¬ 
ternative factorization method. In this analytic method we use the Riccati transformation and 
exponential dichotomies to separate slow from fast dynamics. We employ our factorization 
procedure to study the spectra associated with spatially periodic pulse solutions to a general 
class of multi-component, singularly perturbed reaction-diffusion equations. Eventually, we 
obtain expressions of the slow and fast reduced Evans functions, which describe the spectrum 
in the singular limit. The spectral stability of localized periodic patterns has so far only been 
investigated in specific models such as the Gierer-Meinhardt equations. Our spectral analy¬ 
sis significantly extends and formalizes these existing results. Moreover, it leads to explicit 
instability criteria. 

1 Introduction 

Localized patterns arise frequently in evolution problems with a strong spatial scale separation, 
which naturally leads to the question of their dynamic stability. Several methods have been devel¬ 
oped to study the often decisive spectral stability, especially for the paradigmatic class of parabolic 
semi-linear reaction-diffusion systems on the line of the form, 



( 1 . 1 ) 


where 0 < e 1 is asymptotically small and D\ i are (strictly) positive diagonal matrices. The 
majority of these methods are built on the complex analytic Evans function £ e (T), which vanishes 
precisely on the spectrum and thus is a tool to locate the critical spectrum - see SI1I2Q1 and Remark 
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The singular nature of the equations and the patterns under consideration can reduce the com¬ 
plexity of finding the roots of the Evans function: the slow-fast structure in the linear stability 
problem induces a factorization of the Evans function into a slow and a fast component, 

£ e (A) = £ SiB (A) ■ £ f JA), (1.2) 

as was first observed by Alexander, Gardner and Jones CD in the context of traveling pulses in the 
FitzHugh-Nagumo equations. While operating in a geometric framework, they showed that the 
unstable bundle formed from the projectivized linear stability equations splits into a Whitney sum 
of two line bundles, which are governed by slow and fast singular limit problems. This splitting 
can be directly linked to the factorization < | 1 ,2| > of the Evans function in a slow and a fast com¬ 
ponent. However, although the geometric arguments behind the decomposition in HI are very 
general, they need to be based on an analytical result that indeed controls the relevant subbundles, 
which has been established in lITl l29il in the context of the FitzHugh-Nagumo equations. 

In subsequent work, Gardner and Jones li22ll validated the geometric argument of HI and thus 
the factorization ( |1.2| ), in the context of traveling fronts in a predator-prey model. They tracked 
the fast subbundle analytically via the so-called ‘elephant trunk lemma’. Then, using the control 
on the fast subbundle, they approximated the slow subbundle. Further technical adaptations of 
the elephant trunk le mm a to spectral stability problems associated with localized homoclinic or 
heteroclinic structures in the Gray-Scott and Fabry-Perot model have been developed in lfT3l and 
fl44l . respectively. Nowadays, it is widely accepted that the elephant trunk procedure can be mim¬ 
icked - or better: adapted - for a large class of singularly perturbed systems. However, for every 
application one should in principle go through one of the extensive proofs developed in the setting 
of the aforementioned specific systems to check whether technicalities still hold true. 

By tracking the fast subsystem/subbundle through the elephant trunk procedure, it is possible 
to derive an explicit complex analytic fast reduced Evans function £fo(A) whose zeros approx¬ 
imate those of £/>(/!) - see also Ifl2ll . Moreover, an explicit, but meromorphic, slow reduced 
Evans function £ St o(A) can be obtained via the so-called ‘NLEP (= NonLocal Eigenvalue Prob¬ 
lem) approach’, which was developed in |[T3ll in the context of stability of homoclinic pulses in the 
Gray-Scott model. Thus, a combination of the elephant trunk procedure and the NLEP approach 
yields a reduced Evans function, 

£ 0 (A)=£M-£ ffi (A), (1.3) 

whose zeros approximate those of £ S (A) and whose factors £ Sj o(d) and £/,o(T) can be derived 
explicitly as e —> 0 through properly scaled slow and fast singular limit problems. The NLEP 
approach and thus the validation of the decomposition ( |L2| > and its explicit reduction ( |1 .3[ ), was 
further developed in the context of localized homoclinic or heteroclinic pulses or fronts in certain 
classes of 2- or 3-component singularly perturbed reaction-diffusion equations in lfT2l fT4 , 18:1 55] 
l58l . It should be remarked that in neither of these papers the elephant trunk procedure is worked 
out in full analytical detail. 


This leads us to the main goal of this paper. We establish the validity of both the decomposition 
( |L2| ) of the Evans function £,, : (A) and its singular limit structure (1.3) for a large class of nonlin¬ 
earities H(u, v, s), G(u, v, e) in system <o and general dimensions n,m > 1. Hence, we provide 
a generalized analytic alternative to both the geometric elephant trunk and NLEP procedures. The 
method presented here is based on the Riccati transformation QUO. This transformation, which 
satisfies a matrix Riccati equation, diagonalizes the linear stability problem and thus explicitly 
separates fast from slow dynamics. This separation yields the factorization of the Evans function 
(|L2|) and provides a framework for the passage to the singular limit (| 1 .3[). 
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The factorization procedure can easily be outlined in a way that does neither depend on the specific 
structure of the system nor on the specific patterns under consideration. The first step is to write 
the spectral problem associated with the linearization of <|1.1[> about a singular, localized pattern, 


£s 


u 

V 


f-hduu 

£ 2 D 2 d xx v 


- S e (x) 


u 

V 



as a 2(n + m)-dimensional slow-fast system in block matrix form, 


I d x if \ _ I yfsJlnAx, X) \[s$h2,s(x, X) \/ if \ 

\ d x ft ) [ W 2 i , e (x,T) Xl 2 2 , s (x, A) )\ ft )’ 


(1.4) 


in the rescaled spatial variable x = s ] x. Subsequently, one establishes that the 2/i-dimcnsional 
fast singular limit problem, 


d x ft - Xh2,o(x, A)ft, 


(1.5) 


in which 3\. 2 2 $(x , A) represents the singular limit of J\ 22 ,s(x, 4), has exponential dichotomies on 
both half lines due to the localized nature of the solutions under consideration. These dichotomies 
can be pasted together to an exponential dichotomy on the whole real line, as long as A is not 
an eigenvalue of O- In the third step, roughness techniques are used to carry the exponential 
dichotomy of < | 1 .5| > on R to the perturbed problem, 

d x ft = &22 , £ (x, A) ft, (1.6) 


for 0 < e <sc 1 and A away from the eigenvalues of ( |1.5[ ). This exponential dichotomy on R 
of (JTT6J) allows us to successfully diagonalize the linear stability problem ( | 1 ,4| ) with the Riccati 
transformation yielding a factorization (1.2 1 of the Evans function in a ‘slow’ factor £ StS (A) and a 
‘fast’ factor £/; £ (4). In the last step, we approximate the two blocks, in which ( |1.4[ ) diagonalizes, 
by their singular limits. Once we have understood these singular limit problems, we are able to 
explicitly obtain leading order expressions £. s ,o(d) and £/,o(4). As a consequence, the roots of 
the Evans function can be approximated by the roots of the reduced Evans function £o(4) given 
in 0 Since £ v ,o(d) is meromorphic and £/,o(4) is analytic, one determines the location of the 
spectrum by calculating the roots and poles of £ ij o(^) and £/,o (A)- 


In summary, the complex task of finding the roots of the Evans function £ S (A) associated with the 
2(n+n?)-dimensional system 0 reduces to determining the roots and poles of £ s j)(A) andfi/ oEl), 
which can be derived from lower-dimensional singular limit problems. More precisely, £/,o(A) is 
an Evans function associated with the 2/i-dimensional fast singular limit problem ( |1.5[ ). In contrast, 
£ Vj 0 is an explicit expression in terms of a particular solution to an inhomogeneous version of (1.51 
and the leading order evolution of the 2m-dimensional slow subsystem d x <p = xJeJI n, e (x, A)ip. We 
emphasize that some of the analytic techniques used in the factorization process can be linked to 
the geometric concepts developed in |[Q |22l in the setting of the elephant trunk procedure and in 
fT3l IT211 for the NLEP approach, which is discussed at the end this paper. 


We employ our factorization method to determine the critical spectra associated with stationary, 
spatially periodic pulse solutions to the general class of systems o>- Recall that the proof of 
the elephant trunk lemma has only been worked out in full analytic detail for some specific 2- 
component systems lH3l [T9l 12211441 . Hence, the choice for a large class of multi-component 
reaction-diffusion systems illustrates the general setting to which our method applies. In addition, 
the choice for periodic patterns is motivated by the fact that for certain nonlinearities H(u, v, e), 
there is no simple modification of the elephant trunk procedure that works in this setting. We 
elaborate on the latter claim. 
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It is a general principle that pattern solutions to singularly perturbed models as <o» can be 
‘built’ from exponentially localized pulses (or fronts) in the fast component(s), but allow for non- 
localized behavior of the slow components. In other words, the solutions exhibit semi-strong 
interactions lfT51l (of second order Ii43il ). Most spectral analyses, which make use of the elephant 
trunk and NLEP procedures, are of ‘slowly linear’ nature, in the sense that the dynamics of the 
slow component in between localized fast pulses/fronts are driven by linear equations - see Re- 
mark | 1.1 1 This slow (non)linearity plays a crucial role in the analysis of the Evans function and its 
decomposition and reduction. In fact, it is essential for an application of the elephant trunk proce¬ 
dure to periodic patterns that the matrix in the linear stability problem ( |1.4[ ) is to leading order of 
constant coefficient type near the boundaries of the spatial domain (determined by the periodicity 
of the pattern) - as is the case in lf53l . In slowly nonlinear systems, the matrix ( |1.4[ ) explicitly 
varies in x over the entire domain, thus obstructing an application of the elephant trunk lemma. 

Before applying our factorization method to study the spectra of these periodic pulse solutions, 
we consider their existence. An essential observation is that stationary solutions to o> satisfy 
an ordinary differential equation that admits a reversible symmetry x i—> -x. With the aid of 
geometric singular perturbation theory we establish the existence of reversible symmetric peri¬ 
odic pulse solutions to the (m + //(-component slowly nonlinear system ( | 1.1 [ ). This by itself is 
a significant extension of similar results in the literature that only consider 2-component slowly 
linear Gierer-Meinhardt type models llT6il . However, we emphasize that our spectral analysis does 
not rely on reversible symmetry arguments. In fact, we do not assume that the underlying spa¬ 
tially periodic patterns are reversible. This makes an extension of our methods to models with 
(symmetry-breaking) convective terms or to traveling wave trains natural. 

Our general results can also be interpreted in more simple cases in which either n - 1, m - 1, 
or both n = m = 1. In the latter case, we directly recover the expressions obtained in l53l for 
the spectral stability of spatially periodic pulse patterns in the Gierer-Meinhardt equation. The 
outcome of our spectral analysis shows that the n = m = 1 Gierer-Meinhardt setting represents 
a very special case. This restriction hides the underlying general structure of and Sffi(A) 

in terms of the singular limit problems as obtained here. On the other hand, the restriction of 
<o> to a more general slowly nonlinear 2-component model as in Il8l yields a (relatively) simple 
instability criterion in terms of the signs of a number of explicit integral expressions that can be 
computed with only an asymptotic approximation of the underlying pattern as input. Thereby, we 
extend a similar result of il8l on homoclinic pulses to spatially periodic patterns. We also refer 
to an companion paper in preparation ifTOl . in which the nature of the mechanisms driving the 
destabilization of spatially periodic patterns in general slowly nonlinear systems ( |1.1[ ) - but with 
n - m - 1 - are studied in full analytical (and computational) detail. Thereby, it provides insight 
in the generic nature of the ‘Hopf and belly dances’ discovered in ifTTll in the context of the slowly 
linear Gray-Scott and Gierer-Meinhardt models. 

This paper is organized as follows. In Section [2] we introduce the class of reaction-diffusion 
systems under consideration. Moreover, we elaborate on the existence of periodic pulse solutions. 
Section[3]contains the main results of our spectral analysis. In Section[4]these results are expanded 
further in the case one of the components is scalar. The construction of the Riccati transformation 
is performed in Section[5] The actual spectral analysis via the analytic factorization method is pre¬ 
sented in Section[6] Section [7] contains some concluding remarks and future research possibilities. 
In Appendix [A] one can find the proof of the existence result stated in Section [273] Moreover, in 
Appendix [B] we treat the prerequisites needed for our spectral analysis. Appendix [C] contains the 
proofs of some technical, but not fundamentally difficult, results in this paper. 
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Remark 1.1 (Slow nonlinearity). As mentioned earlier, most spectral analyses, using the elephant 
trunk and NLEP procedures, are done in models of ‘slowly linear’ nature, which include the 
classical Gray-Scott and Gierer-Meinhardt models. In recent work Ifl8l l58ll an NLEP approach 
has been developed for the spectral analysis of homoclinic pulses to a general class of singularly 
perturbed slowly nonlinear 2-component reaction-diffusion systems. Earlier, spectral stability of 
fronts was studied in a specific slowly nonlinear model in Ifl4ll . Although the models in these 
works are slowly nonlinear, the elephant trunk procedure is still applicable, because eventually 
the dynamics of the slow component becomes linear due to homoclinic or heteroclinic nature of 
the patterns. However, we emphasize that there are significant adaptations necessary in the NLEP 
procedure for an extension of the slowly linear to the general slowly nonlinear case - see ifUfl . ■ 


Remark 1.2 (Spectral stability). It should be noted that only studying the roots of Sffi(A) and 
Sv.of/l) is not sufficient to decide upon spectral stability of the spatially periodic patterns, because 
there must be essential spectrum attached to zero: the leading order results of the present work 
do not provide information on the exact position of this ‘small spectrum’ ll53l with respect to the 
imaginary axis. Therefore, an additional study of the fine structure of the spectrum around zero is 
necessary to determine spectral stability - see Section [374] This is subject of work in progress. ■ 


Remark 1.3 (The Evans function and other tools to locate the spectrum). The concept of the 
Evans function as a method to determine the spectrum associated with a localized solution of 
a system of reaction-diffusion equations on the line was introduced in Il20l and was established 
as a general and powerful approach in [1] [221 [29]] • Core aspects of the NLEP approach have 
been developed independently in lf28l l59l . The SLEP (= Singular Limit Eigenvalue Problem) 
method Il36l [37 1 is an alternative method that has been linked to the Evans function approach in 
1271 . In |H, the relation between the Evans function, the NLEP method and the SLEP method 
is discussed. The Evans function approach was originally developed in the context of localized 
homoclinic and/or heteroclinic patterns, it was first generalized to spatially periodic pulse patterns 
in reaction-diffusion equations in fl9ll23ll53ll . ■ 


2 Setting 

2.1 Reaction-diffusion equations with semi-strong interaction 

In this section we introduce the class of systems under consideration in this paper. Take m, n e Z>o 
and consider a general reaction-diffusion system in one space dimension with a scale separation 
in the diffusion lengths ( [ 1.1 1 ). Following lfT8l . we write 

H(u, v, s) = H(u, 0, s) + Hfu, v, s), 

where Hfu, v, e) := H(u, v, e) - H(u, 0, e), so that Eh_ vanishes at v = 0. To sustain stable localized 
patterns in semi-strong interaction (of second order 1431 ) in system <0>’ we allow Hfyu, v, s) to 
scale with e 1 and define 


Hi(u,v) := lim eH->(u,v,e). 

sj.0 


Finally, we write 


H(u,v,s) = H[(u,v,s) + e 1 H2(u,v), 


(2.1) 


with Hi(u,v,s) := H(u,0,e) + [TLfffiv, s) - e 1 //?(«, v)]. By construction H 2 (u,v) vanishes at 
v = 0. Moreover, we assume that Hfu, v, s) and G(u, v, s) are smooth functions of e at s = 0. Note 
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that we allow for the possibility that H 2 (u, v) = 0 in the upcoming analysis. We emphasize that, 
when both H 2 (u,v) = 0 and n = 1, all patterns are unstable - see Remark 3.13 This confirms the 
scalings used for classical systems as the Gray-Scott and Gierer-Meinhardt models lfT3l IT2l 12811591 
- see also lfT8l . For the benefit of our spectral analysis, we need one extra condition on G. That 
is, G vanishes at v = 0. We postpone the discussion of this extra condition to Remark 2.6 In 
summary, the model class we consider is of the form 


d t u = D\d xx u — H\(u,v,s) — s l H 2 (u,v) 
d t v = s 2 D 2 d xx v - G(u, v, s) 

or, in the ‘small’ spatial scale x = s ~ 1 x, 

e 2 d t u = D | d xx u - e 2 H\ (u, v, s) - eH 2 (u, v) 


u € R. m , v e R", 


( 2 . 2 ) 


d t v = D 2 d xx v - G(u, v, s ) 
in which we will usually work. The aforementioned conditions read: 


u e R , v e 


(2.3) 


(SI) Conditions on the interaction terms 

There exists open, connected sets U c R m , V c R" and I c R with 0 e V and 0 € / such 
that H\,G and H 2 are C 2 on their domains U xV x I and U x V, respectively. Moreover, we 
have H 2 (u, 0) = 0 and G(u, 0, e) = 0 for all u e U and eel. 


Remark 2.1 (The quantities D\ 2 in the cases m = 1 and n = 1). If we have n = 1, we can without 
loss of generality assume D 2 = 1 in (2.2) by rescaling the spatial variable x. Similarly, in the case 
in — 1, we can without loss of generality assume D\ = 1 by rescaling the small parameter e. m 


2.2 Stationary, spatially periodic pulse solutions 


In this paper we are interested in the spectra associated with stationary, spatially periodic pulse 
solutions satisfying the reaction-diffusion system (2.3). Stationary solutions to (2.3) satisfy the 
singularly perturbed ordinary differential equation, 


D\d x ti 

d x p 

D 2 d x v 

d x q 


sp 

eH\(u , v, s) + H 2 (u, v) 

q 

G(u, V, E ) 


ue U,p er,ve V,q e R”. 


If we take s = 0 in (2.4 1 , the dynamics is given by the so-callcd/r/.v? reduced system, 


d x u 

d x p 

D 2 d x v 

d x q 


0 

H 2 (u, v) 

q 

G(u, v, 0) 


ue U,pe P,ve 


(2.4) 


(2.5) 


We observe that W4 = ((«, p, 0,0) : u € U,p e IF"' j consists entirely of equilibria of ( |2.5| ) by 
assumption|(Sl)| We require M to be normally hyperbolic. 


(S2) Normal hyperbolicity 

For each u e U the real part Rc(Q(u)) = 
definite. 


5 (Q(u) + Q(u)*) of Q(u) d v G(u, 0,0) is positive 
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When e > 0, the manifold At consists no longer of equilibria, but remains invariant. The flow 
restricted to At is to leading order governed by the so-called slow reduced system , 


D\dxU 

dxP 


P 

H\{u, 0,0) 


u e U. p e R m . 


( 2 . 6 ) 


It is well-known that the dynamics around such a normally hyperbolic manifold At is captured 
by FenicheTs geometric singular perturbation theory li2Tl [30 1. Fenichel’s theory can reduce the 
complexity of finding periodic orbits in the following way. Suppose we have obtained a so-called 
singular periodic orbit by piecing together orbit segments of the fast and slow reduced systems 
in such a way that they form a closed loop. Although this singular orbit is not a solution to the 
full system, one can prove with the aid of Fenichel’s theory that an actual periodic orbit lies in the 
vicinity of the singular one, provided s > 0 is sufficiently small. This construction is performed in 
Section [2731 


In this paper we are interested in periodic solutions to ( |2.4| ) that are close to singular periodic 
orbits, which consists of two orbit segments, to wit: a pulse satisfying the fast reduced system 
(2.5 i and a segment on the invariant manifold At, satisfying the slow reduced system (2.6). We 
make the hypotheses on the solutions precise in the following assumption. 


(S3) Existence of a periodic pulse solution 

For e > 0 sufficiently small, there exists a 2L e -periodic solution <p V E : R —> U x R m x V x R' ! 
to system (|2.4|) with the following three properties. 


1. Order of the period 

It holds L e = e~ 1 L £ and L s —» Lo > 0 as s —> 0. 

2. Existence of a pulse 

There exists uq e U, po £ R m and a solution (vh(x), qAx)) to system. 


D 2 d x v 

d x q 


q 

G(uq, v, 0) 


ve V,q e R' ! , 


which is homoclinic to (0,0), such that 


H0p,e(O) - (n 0 , Po, Vh(0), #h(0))|| = O(e). 


(2.7) 


( 2 . 8 ) 


3. Exponential decay to a solution in At 

There exists L S ,/3 S > 0 and a solution (w s (x), pfx)) to (2.6) such that it holds for 

x € [0,2L £ ] 


d(4>p,s(x). At) < min(x ’ 2Le - Jc| , (2.9) 

||0p,e(x) - (uJex), Ps(ex), 0,0)|| < L. v max {e, e~^ s mm A,2L e -x\\ _ (2.10) 


Remark 2.2 (Reversible symmetry). Note that (2.4) is irreversible, where R : R 2 ( m+ ") R 20 " + ") 

is the reflection in the space p = q = 0. Indeed, if 0 is a solution to ( |2.4[ ), then x Rcp(-x) is 

1 2m _> r 2»> and R f 


also a solution to (2.4 1 . Similarly, denote by R s : 


V- 


»2 n 


s2 n 


the reflections 


in p = 0 and q = 0. respectively. Systems (2.6) and (2.71 are R s - and i?j-reversible, respectively. 


Reversible symmetry arguments significantly simplify our existence analysis in Section 2.3 As a 


consequence, the constructed periodic orbits are symmetric about the space p = q = 0. However, 
reversible symmetry arguments do not essentially simplify our spectral analysis. In fact, these 
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arguments could obscure an extension of our results to models with convective terms, which do 
not obey the reversible symmetry. Therefore, we do not impose a reversible symmetry assumption 
on the solutions (p VE in assumption |(S3)| ■ 


Remark 2.3 ( Embedding in an E-independent compaction). An immediate and important conse¬ 
quence of (S3) -3 is that there exists a compact set 7C c U x R m x V x R", independent of e, such 
that (pyij:(x) 6 r K for all x 6 R. Note that we can choose 'K in such a way that its projection 'Ku c U 
on the ^-component contains uq. This will be convenient for later purposes. ■ 


Remark 2.4 (Exponential decay). By assumption (S2) (0,0) is a hyperbolic saddle in system 


(2.7 1 . Hence, the homoclinic (vh,r/h) is exponentially localized, i.e. there exists constants K,p > 0 
such that ||(vh(x), <7h(-T))ll ^ Ke~^ for x e R. ■ 


Remark 2.5 (Extending estimate (2.8 1). Consider the solution <!>b(x) to system (2.5 I given by 
(,K{x ) = |mo, Po + H 2 (u 0 , v h (z))dz, v h (x), <?h(A>j • 


System ( |2.4[ ) can be written as d x if/ - f(ip, e), where / is locally Lipschitz by |(S1)| Denote by Lf 
the Lipschitz constant of / on 7C x [0, eo] for some 0 < eo € 1. We infer via Gronwall’s inequality 
and assumption |(S3)[- 2 that there exists a constant K > 0, independent of e, such that 


ll0p, £ W - 0h(*)ll < e(K + Lf\x\)e L f ]x] , 


( 2 . 11 ) 


for jc 6 R. First, we employ ( |2. 11 [ ) on J e := [x\ log(e), —x\ log(e)] with x\ > 0 e-independent. 
Second, we use Remark |2.4| and (S3)} 3 to derive estimates of d(cp v ^ E (x), M) and dUp\dx), M) on 
[- L £ ,L e \ \ J E . Putting these items together yields /?/,p > 0 such that 

UpAx) - 0hMII = 0(eP f), 

for x e I E := [~£~ p , s ~ p ]. We will refer to I E as the pulse region. ■ 


Remark 2.6 (The condition that G vanishes at v = 0). As mentioned in Section |2. 1 our model 
(2.2 1 is a general reaction-diffusion system that allows for semi-strong interaction fl2.1| ), with the 
extra condition that the term G vanishes at v = 0. For general G, consider a 2/z-dimensional com¬ 
pact submanifold Alo of {(u, p, v, 0): G(u, v, 0) = 0} c U X R m x V x R”. By Fenichel’s theory Ato 
perturbs, for e > 0 sufficiently small, to a locally invariant manifold Al e to (2.4). This manifold 
Al £ is diffeomorphic to Alo and lies at Hausdorff distance (9(e) from Alo. When Alo can be given 
as a graph over (; u , p) € U x R m , the same holds for Al £ . Thus, in that case one can change coor¬ 
dinates in \2A\ relative to Al e and we obtain Al £ = Alo c Al. Therefore, on the existence level, 
the condition that G vanishes at v = 0, corresponds to an a priori coordinate change in (|2.4 1 . 


However, on the stability level, one introduces more than additional technical difficulties when 
G does not vanish at v = 0. Indeed, without relative coordinates, we do not achieve estimate ( |2.9[ ), 
which is essential in our spectral stability analysis. However, applying the coordinate change to 
equation ( |2.2[ ) changes its structure fundamentally. In the new coordinates ( |2.2[ ) is not even of 
reaction-diffusion type. Hence, we expect that the spectral analysis differs essentially, when G 
does not vanish at v = 0. This is an interesting subject of future research, especially since it in¬ 
cludes the possibility of localized patterns with oscillating tails (6j[T8J, but is outside the scope of 
this paper. ■ 























2.3 Existence of stationary, spatially periodic pulse solutions 


To motivate assumption |(S3)[ we elaborate on the existence of reversible periodic pulse solutions 
to system \2A) satisfying assertions (S3) 1,2,3. For a general class of slow-fast systems the fol¬ 
lowing result is established in ll52ll . Suppose a singular periodic orbit is given, obtained by piecing 
together orbit segments of the fast and slow reduced systems in such a way that they form a closed 
loop. Then, under certain conditions, an actual periodic orbit lies close to the singular one, pro¬ 
vided e > 0 is sufficiently small. 


An essential condition for the result in lf52l is that the slow components are constant in the fast re¬ 


duced system. In our case the slow /;-component is non-constant along orbits in ( |2.5[ ). Therefore, 
the result in ll52l is not applicable. However, we will obtain a similar result. Consider a singular 
periodic orbit (Figure [Tb|) obtained by gluing a pulse solution to (2.51 to a solution on At governed 


by (2.6 1 . For small e > 0, we will show that there exists a periodic orbit in (2.4 1 in the vicinity of 
this singular periodic orbit. 


Concerning the existence of the singular periodic orbit, we need some assumptions. We start 
with the assumption on the fast reduced system (|2.5[). 


(El) Existence of a pulse solution to the fast reduced system 

There exists uq e U such that (2.7) has a solution fh(x,uo) = {\\{x,uo),q^{x,uf)) homo- 
clinic to 0. The stable manifold (0) to the hyperbolic saddle 0 in system (2.7 1 intersects 


the space ker(7 - Rf) transversely in the point fh(0, uq), where Rf is as in Remark 2.2 


Remark 2.7 ( Elementary homoclinic orbits). In the terminology of lf56l homoclinics that cross 
a transverse intersection of VF,] (] (0) and ker(7 - Rf) are called elementary. In particular, any non¬ 
degenerate homoclinic solution is elementary by ll56l Lemma 4]. ■ 


Remark 2.8 ( Extension of assumption \(E1]\ . Since transverse intersection are robust under per¬ 
turbations, assumption |(E1)| implies the existence of an open neighborhood Uh c U of uq such 
that for every u e Uh there exists a homoclinic solution fh(x, u ) to 


D 2 d x v 

d x q 


<7 

G(ii, v, 0) 


v € V,q e R", 


crossing (0) iti ker(7 - Rf) = j///h(0, «)!■ 


Remark 2.9 ( Take-off and touch-down manifolds). The homoclinics Jh(x, u) for u e Uh yield 
solutions, 




(p h(x, u) := w, I 77 2 (m, V h (z, u))dz, v h (x, u), q h (x, u) , 


to (2.51, which are homoclinic to At The limits lim AW±0o 0h(r > Ll ) give rise to the so-called take-off 
and touch-down manifolds. For that reason, we define the mapping ff : Uh —» R' H by 


-f 


J(u) = H 2 (u,vh(z,u))dz. 


The m-dimensional graphs T± := ((«, ±ff(u)) : u € 77h} on Al are the take-off and touch-down 
manifolds. Indeed, it holds lim r ^ ±00 (pjx, u) = (u. ±ff{u), 0,0). ■ 


9 























The manifolds T ± of Remark |2.9| allow us to piece the pulse solutions (t>h to solutions that lie in 
M in order to obtain a singular periodic orbit (Figure |Th| ). Therefore, we shift our attention to 
the slow reduced system (2.6). Recall that system \2.6\ is 7? s -reversible by Remark |2.2| It holds 
R S VT+ 1 = 71. Hence, to establish a connection between the take-off and touch-down manifolds 
71, it is sufficient to find a solution that starts in ker(7 - R s ) and crosses the take-off manifold 71 
at some point. This is the content of our next assumption. 


(E2) Existence of connecting orbit in slow reduced system 


Let R s ,T±,ff be as in Remarks |2.2| and |2.9[ There exists a solution ffx) = (u s (x), p s (x)) to 
system (2.6) such that i// s (0) e ker(7 - R s ) and if/fLf) e 71 for some Lq > 0. Moreover, let 
cD s (x,y) be the evolution of the associated variational equation, 


d.x<p - 


0 


D 


,-i 


duHfufxfO'O) 0 


V- 


Denote for /, j e {1,..., m} by hj the coordinates of H\ (u s (0), 0,0) and by Ajj the (m x /;?)- 
submatrix of 


O s (0,L 0 ) 


-duff («o) 


containing rows {/, m + 1,..., 2 in) \ {m + j }. There exists ;* € {1,..., m) such that 

m 

2(-l)4jdet(A w )*0. (2.12) 

7=1 


To obtain persistence under small perturbations, one could require that the intersection between 
the take-off manifold 71 and the trajectory is transversal. However, this is impossible for m > 1. 
Therefore, the technical condition ( |2.12[ ) is employed to generate a ‘good’ set of initial conditions 
in ker(7 - R s ). This set becomes under the forward flow of the slow reduced system ( |2.6[ ) an m- 
dimensional manifold, which contains the solution and does intersect 71 transversally (within 
M). 


Remark 2.10 {Assumption \(E2)\ in the case m = 1). In the case m = 1 system (2.6 1 is Hamiltonian. 
Therefore, assumption |(E2)| is satisfied if and only if there exists u\ e U such that the following 
four conditions hold true. 


f 

*Ju\ 


i 


Hi(u,0,0)du = -ffiuff, (u\ 


uo)J(uq) > 0, 
77i(«i,0,0) + 0, 77i(m 0 ,0,0) + J(u 0 )d u J(u 0 ). 


The first three conditions can be related to the existence of a trajectory from (u \, 0) e ker(7 - R s ) 
to («o, ~ff(u o)). The last condition is associated with the transversality of the intersection between 
this trajectory and the take-off curve 71. ■ 


We are now able to state the existence result. The proof uses techniques developed in fl6l . where 
the existence of stationary, spatially periodic pulse solutions in the Gierer-Meinhardt equation 
is considered. However, we emphasize that the framework in lll6l differs from ours due to a 
difference in scaling in the p-component. This is explained in more detail in Remark [3Tj 


10 














(a) Slow-fast dynamics: the v-component exhibits local¬ 
ized pulse behavior and the n-component evolves slowly. 



(b) Orthogonal projection of the sin¬ 
gular periodic orbit on the slow man¬ 
ifold M. 


Figure 1: Periodic pulse solutions in the case n - m- 1. 


We extend and adapt the techniques in Ifl6l to the more general existence problem (2.4). This ex¬ 
istence result, which focuses on periodic solutions in a large class of singularly perturbed systems 
beyond the type of slow-fast systems considered in ll52ll . is to our knowledge new - we are not 
aware of a similar general result in the literature. To avoid too much diversion, we present the 
proof in Appendix [A] The proof has a geometric nature. It relies on Fenichel’s theory, reversible 
symmetry arguments and an appropriate version of the Exchange Lemma. 


Theorem 2.11. Assume [(SI ) (S2) and (El) (E2)\hold true. Then, for s > 0 sufficiently small, 


there exists a 2L s -periodic solution 0 PiS to (2.4 ) satisfying assertions 1-3 in (S3) Moreover, f p _ £ is 


symmetric about the space p = q = 0. It holds f p , £ (x) = Rf p£ (-x) and <p p ^ £ (L £ -x) = Rf p£ (L £ + x) 
for x e R, where R is as in Remark 2.2 


Remark 2.12 (An expression for Lo in th e ca se m = 1). Let m = 1 and ti\ e U as in Remark 2.10 
Using the Hamiltonian nature of system (2.6 1 , we can express Lo, as defined in (S3) T, in terms of 
mq and mi as 


Ln — 


rn o 
Jll\ 


du 


ui sj2F(u) + J(u 0 ) 2 ' 


where F(u) := H[(u,0,0)du for u e U. 


Example 2.13 (Generalized Gierer-Meinhardt Equations with slow nonlinearity). An instance of 
( |2.3| ), which satisfies assumptions (SI) (S2) and (El) (E2)[ is the generalized Gierer-Meinhardt 
equation with slow nonlinearity, 


d t u 

d t v 


sd xx ii - erf(u) + su ai v@ l 
d xx v - v + u ai vP 2 


U E ]R>o> V E R, 


(2.13) 


with parameters a\^ £ R,/3i, 2 £ Z>i and / € C 2 (0, oo). One can take f(u) = pu for some p > 0 to 
obtain the slowly linear Gierer-Meinhardt equations considered in If53l . In the model considered 
in ll58l one takes f(u) = c\u - C 2 U d with c\p > 0 and d > 1. Also, the choice f(u ) = sin(n) 
is possible. Then, one obtains the sinusoidal Gierer-Meinhardt model considered as a guiding 
example in the companion paper iflOl . ■ 
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3 Main results 


Throughout the rest of this paper, assume (S1)||(S3) and regard <p v , e (x) = (n P;£ (x), v P;£ (x)), for 
e > 0 sufficiently small, as a stationary, spatially 2L £ -periodic pulse solution to (2.3) as de¬ 
scribed in (S3) Let 0 Pi£ (x) = (w p , £ (x), v p , £ (x)) be the corresponding solution to the rescaled system 
(2.2). We linearize system (2.2 1 about r/» p £ and obtain the differential operator £ s : C 2 b { R, I 
n+n ) -» C ub { R, R m %given by: 


C ub {l 


with 




D\dxxM 

£ 2 D 2 diiV 


- 8 e (i) 


S e (x) := 


d u Hi((p PtS ,s) + e l d u H 2 (<f > p , e ) d v Hi($ p , s ,s) + s l d v H 2 (0 p , 6 ) 
d u G{(j) p, £ , £) d v G(<f> p , £ , s) 


) c 


(3.1) 


(3.2) 


suppressing the ^-dependence of r/> p £ in the right hand side of ( |3.2| ). Here, CL(R, R m+ ”) denotes 
the Banach space of k times continuously differentiable functions, with derivatives up to order k 
bounded and uniformly continuous. It is endowed with the norm. 


= £ IKWIU- 


i '=0 


Note that by assumption (S3) l, X £ is a periodic differential operator with period L £ . 


In this section we collect the main outcomes of our spectral analysis. We formulate a leading 
order approximation result for the spectrum c t(£ e ). This leads to explicit (spectral) instability 
criteria. As mentioned in Remark [L2] the leading order approximation results are not sufficient to 
decide upon spectral stability. However, we will show that the question of spectral stability can be 
reduced to a study of the fine structure of the spectrum around 0. 


3.1 Formulation of the Evans function 

Before stating the main results, we focus on the structure of the spectrum cr{£ £ ) and introduce the 
analytic Evans function that vanishes precisely on cr(X £ ). 


Regard Xc, £ as the extension of £ E to C ub { R, C m+ ") given by 


£c,e<P = £*[Re(<p)] + LC e [Im(^)|. 


Fundamental in our spectral analysis is that by Il23l Proposition 2.1] a point A e C is in the 
spectrum cr(X £ ) if and only if there exists <p e C 2 h ( R, C m+ ")\!0) such that Xc = d</>. In the small 
spatial scale x = £~ 1 x the eigenvalue problem £c,e^P = can be written as the linear system, 


with periodic coefficient matrices, 


(3.3) 


ftn,e(x>X) '■ = 
&12A x ) : ~ 


D 


,-t 


£ (d u H\ (w Pj£ (x), v p , £ (x), E) + d) + d u H 2 (Up^(x), v p , £ (x)) 0 J ’ 

0 0 \ 

£d v Hi(u PtE (x), v P;£ (x), e) + d v H 2 (u v Jx), v p , e (x)) 0 / ’ 
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and 


^ 21 , six) := 


0 


0 


di,G(u„ s {x),v p, E (x),e) 0 


^22,eU, d) : = 


0 


D 


,-i 


d v G(llnJx), V P)£ (x), s) + A 0 


We will refer to (3.3) as the linear stability problem. 


Remark 3.1 (Scaling of the p-component). We emphasize that there is a difference in scaling in 
the p-component between the linear stability problem (3.31 and the existence problem \2A\. The 


scaling regime in (3.31 brings the system in slow-fast form, which has the advantage of making 
important tools available such as the Riccati transform (Theorem |5.1[ ). In principle, the existence 
problem could also be put in this form by introducing p - £~ l/2 Did x ii instead - as is done in exis¬ 
tence analysis of periodic pulse solutions in the Gierer-Meinhardt equations in Iil 6 l . The equation 
for the 72 -component in the slow reduced system (|2.6|> would read d x p = 0 in that case. This 


makes the construction of the desired singular periodic orbit, performed in Section |2.3| impos¬ 
sible. Therefore, the scaling regime in (2.4 > is the most natural for our existence analysis. In 


lfT6l one avoids setting e = 0 in the existence analysis and makes a distinction between slow and 
‘super-slow’ behavior. ■ 


We are looking for solutions e C„/,(R, C" !+ ") to (3.31. It is well-known by Floquet Theory [|4[ 
Chapter 1] that bounded solutions to ( |3.3[ ) must satisfy tp(-L E ) = ytp{L E ) for some y in the unit 
circle S 1 . This fact gives rise to the following definition, which provides a tool to locate the 
spectrum of 


Definition 3.2. Denote by T E (x, z, A) the evolution operator of system (3.3) with initial condition 
z. The Evans function S E : C x C —> C is given by 

G e (A, y) := detCT e (0, -L e , A) - yT E (0, L s , A)). 

Proposition 3.3. The Evans function has the following properties. 

1. The Evans function is analytic in both A and y. 

2. We have A e cr(E F ) if and only if there exists y e S 1 such that the dispersion relation 
S s (A, y) = 0 is satisfied. 


3. It holds Ej,(A, y) = G E (A, y) for A, y e C. Thus, the spectrum <r(f E ) is closed under complex 
conjugation. 

Proof. The first two properties are the content of Il23l Lemma 2.4] and Il23l Proposition 2.1], 


respectively. Since (3.31 is a real valued problem for A € R., the third property follows by the 
reflection principle. □ 

Proposition |3.3| shows that the spectrum of X e is the union of the discrete zero sets N(G e (-,y)), 
where y ranges over S 1 . Hence the spectrum is an at most countable union of curves, each 

of which is parameterized over S [ . This gives rise to the definition of y-eigenvalues, which is due 
to Gardner Il23l . 

Definition 3.4. Let y € S 1 . A zero A of the Evans function £ e (-, y) is called a y-eigenvalue of f E . 


Remark 3.5 (Reversible symmetry). Recall that system (2.4) is A’-rcversihlc by Remark |T2| Sup¬ 
pose obeys this reversible symmetry, which is the case for the solutions obtained in existence 
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result Theorem 2.11 Then, the linear stability problem (3.31 is irreversible at 0, i.e. it holds 
RT E (x,y,A)R = T E (-x,-y,A ) for x,y e R. One readily deduces G E (A,y) = G g (A, y - ')y 2(m+ «) for 
A 6 C and y € C\{0). So, if A e C is a y-eigenvalue for some yeS 1 , then A is also a y-eigenvalue, 
since we have y = y~ l for y e S l . Therefore, the image of S 1 covers each curve of spectrum 
twice, which yields degenerate spectrum (Figure [2]). ■ 


3.2 Approximation of the spectrum 

The main outcome of our spectral analysis is an explicit reduced Evans function Gq(A, y), whose 
zeros, for y restricted to .S' 1 , approximate the zeros of the Evans function G E (A, y), provided that 
s > 0 is sufficiently small. Therefore, the reduced Evans function Go proves to be a powerful 
tool to study the shape of the spectrum cr(X £ ) and decide upon spectral (in)stability questions. 
We emphasize that the approximation of G e by £o is only valid on certain half planes Ca, which 
contain the critical spectrum of _£ e . 


Notation 3.6. For every A < 0 we denote by Ca the open half plane, 

C A := fieC: Re(T) > A}. 


The slow-fast structure of the linear stability problem (3.31 is reflected by the fact that £o can be 
written as a product 


G 0 (A, y) = (~y) n Gy 0 (A)G sfi (A, y). (3.4) 

Here, the analytic map Gy o'- Ca —» C is called the fast reduced Evans functions. It is associated 
with the homogeneous fast limit problem, 

* « c2 "' =- ( a,,G(»„.v h °w.o) + . 1' ) ■ <*» 

System ( |3.5[ ) corresponds to the homogeneous variant of the fast layer problem in the singular 
limit of (|3.3[). We establish the existence of the fast reduced Evans function. 


Proposition 3.7. There exists A < 0 and an analytic map £/o: Ca —> C, which has a zero if 
and only if problem (|3.5|) admits a non-trivial bounded solution. In particular, A € Ca is a simple 


zero ofGyo if and only if there exists a unique non-trivial bounded solution to (3.5 1 up to scalar 
multiples. 


The slow reduced Evans function £ Vi o: [Ca \ £^(0)] x C —> C is determined by two problems. 
The first is the inhomogeneous fast limit problem, 


d x X = ^22,o0t A)X + <^2i,o(x), X e Mat2„ X 2m(Q. ^ 21 ,oW 


0 0 

d u G(u 0 , v h (x), 0) 0 


(3.6) 


The matrix system (3.6 1 can be seen as a family of inhomogeneous fast layer problems in the 
singular limit of (|3.3 1 . The second is the slow limit problem, 


I D\d x u = p 

\ djtp = (d u Hi(u s (x),0,0) + A)u ’ U,P 


(3.7) 


which correspond to the slow subsystem of (3.3). 


14 















Definition 3.8. The slow reduced Evans function £ v .o: [Ca \ £^(0)1 x 


is defined by 


£,,o(d, r) = det (T(d)r s (2L 0 ,0, A) - yl), 

where Lq is given by assumption |(S3)f l, T s (x,y,A) is the evolution of the slow limit problem (3.7 1 
and Y(/l) is given by 


T(d) = 


/ 0 

GW / 


I 


(3.8) 


GW = [d u H 2 (u o, i’h(x)) + d v H 2 (uo, Vh(x))'Vm(x, A)\dx, 


where r Vi n (x, A) denotes the upper-left (n x /7?)-block of the unique (2 n x 2/n)-matrix solution 
X m (x, /l) to the inhomogeneous fast limit problem 


Wc collect some properties of the slow reduced Evans functions £ lSj q. 


Proposition 3.9. The slow reduced Evans function £ v o is analytic in both A and y. Moreover, 
£ v o(', y) is meromorphic on Ca in such a way that the reduced Evans function £o(-,y) is analytic 
on Ca for each y <E C. 


Concerning zeros of £o, consider a simply closed curve Y c Ca that avoids the roots of S s p and 

fi/.o 


M) := 


Af(£/,o) u [J Af(£ J>0 (-,y)) 

yeS 1 


Let y € 5 1 . Denote by AY, AY. y and Pr y the number of zeros of £yo, £ iV ,o(-, y) and poles of 
£.v.o(', y) interior to T counted with multiplicity, respectively. Given the analyticity and product 
form ( |3.4| ) of the reduced Evans function £o(-,y), its number of roots interior to Y counted with 
multiplicity is given by 


N r + Nr,y-Pr,y. (3.9) 

Our first approximation result confirms this and its persistence for £ e with small e > 0. 


Theorem 3.10. Take a simple closed curve Y contained in Ca \Nq. For e > 0 sufficiently small, 
the number ofy-eigenvalues (including multiplicity) of E, : interior to Y equals 


Theorem 3.10 shows that the set, 


o-o ■= [J Af(£o(-,y)X 

yeS 1 

can be seen as the spectrum in the singular limit of the operator E, : in Ca- Indeed, choose contours 
close enough to and disjoint from the connected components of cr 0 , with, say, Hausdorff distance 
5. This results in an > 0 such that, if s e (0, e^), then the spectrum cr(X £ ) n Ca is contained in 
a d-neighborhood of cr 0 . In that sense the reduced Evans function £o determines the shape of the 
spectrum cr(X e ) up to connected components of the set cr 0 . However, this approximation could 
still be too rough, since it does not concern pointwise distances of spectra. For instance, when a 
connected component of cr 0 transversally intersects the imaginary axis, one can not decide with 


Theorem 3.10 upon spectral (in)stability. In that sense, our second approximation result is more 
refined. 
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U 7 eSn£>( - >7) HO) 

(a) Spectrum in the singular limit with disjoint 
contours I | and Y 2 . 

Figure 2: approximation of the spectrum cr(X e ) 


<■ J (£e ) 

(b) The true spectmm remains in the interior of 
Ti and Y 2 for s > 0 sufficiently small. 

in the case of reversible symmetry (Remark|3.5[). 


Theorem 3.11. Fix y e S 1 . There exists 8 y > 0, such that for every 0 < 8 < 8 y , provided e > 0 is 
sufficiently small (with bound dependent on 8), the y-eigenvalues of f, : in C,\ are contained in a 
8-neighborhood of N(Bq(-, y)). Moreover, for every A y e N(So(-,y)) the number of y-eigenvalues 
(including multiplicity) inside B(A y , 8) is equal to the multiplicity ord| (£o(-, y)) of A y as a root of 

fio(-.r)- 


Corollary 3.12. The solution 0 p £ is spectrally unstable for e > 0 sufficiently small, if one of the 
following is true: 

1. There exists y,, e .S' 1 and A,, e C with Rc(/Lj > 0 satisfying <5 Vj o(/l®, y<>) = 0; 

2. There exists y,, e .S' 1 and a zero A 0 of Bp 0 with Re(/t 0 ) > 0 such that ord,{„(£ v .o(£ y,,)) + 

ord^(£/ j0 ) > 0. 


When a connected component of cr 0 transversally intersects the imaginary axis, Corollary 3.12 
shows that the periodic pulse solution is, for e > 0 sufficiently small, spectrally unstable. 
Indeed, there must be in that case a point, which has positive real part, in one of the discrete sets 
N(B s fi(-,y)) for some y e S 1 . The fact that this point can be isolated, is the strength of Theorem 


3.11 The trade-off is that the bound on the small parameter e is y-depcndcnt in Theorem 3.11 


whereas the bound on e is uniform in y in Theorem 3.10 


The proofs of our main approximation results Theorems 3.10 and 3.11 can be found in Section[6j 


Remark 3.13 ( Weak coupling). As mentioned in section |2.1[ weak coupling H 2 (u, v) = 0 is al¬ 
lowed in our spectral analysis. In that case the integral term Q(A) in ( |3.8[ ) is identically 0, which 
implies that B s p is only determined by the slow limit problem < ]3 .7 [ ). Therefore, £ 4j o is analytic on 
Ca and zeros of £y o can not be canceled by poles of £ i; o- We conclude that the spectral stability 
problem fully splits into slow and fast subproblems with no interaction between them. As a con¬ 
sequence, zeros of Bp 0 of positive real part yield spectral instability immediately by Proposition 


3.12 -2. In particular, in the case n = 1 the fast reduced Evans function Bpo always has a zero 
in the right half plane, as we will show in Proposition |4.1| Thus, all patterns are spectrally 


unstable in the case H 2 (u,v) = 0 and n - 1. This motivates the scaling in (|2. 1 [). 
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Remark 3.14 ( Translation invariance and spectral stability). If there is a zero A y of £o(-, y) sat¬ 
isfying Re(/l y ) > 0 for some y € S 1 , we know that the periodic pulse solution 0 p £ is spectrally 
unstable by Corollary |3.12[ However, to determine whether <p PiE is spectrally stable, it is not suffi¬ 
cient to know the zeros of £o. This is due to translational invariance, i.e. the derivative of is 
always an eigenfunction of £ £ . Hence, 0 will always be a zero of £,.(■, 1). So, to decide whether 
is spectrally stable, a study of the so-called small spectrum around 0 together with the fact 
that for all y € S 1 there are no zeros A y + 0 of £ 0 (-, y) satisfying Re(d y ) > 0 affirms us that <p p>E is 
spectrally stable. We refer to Corollary |3.17| for a more rigorous statement. ■ 

Remark 3.15 ( Bloch-wave transform). It is also possible to consider the spectrum of the operator 
£ £ in the Sobolev space of y-twisted periodic functions for fixed y e .S' 1 . The spectrum is in 
that case given by the discrete set A/ r (£ e (-,y)) of y-eigenvalues of The decomposition of 
the spectrum of £ £ in parts AT£ £ (-, y)), where y ranges over S 1 , corresponds to the Bloch-wave 
transformation of the operator £ £ . ■ 


3.3 Zero-pole cancelation 

Take y e S 1 . By Proposition 


3.9 


a pole of the slow reduced Evans function £ v .o(c 7) can only 
occur at a zero of the fast reduced Evans function Gf o, because the reduced Evans function 
£o(/l, y) = (-y)"£/,o(/l)£ s ,o(d,y) is analytic. In particular, by formula (3.91 in Theorem 3.10 it 


can occur that £yo has a zero A e C a with positive real part, but the periodic pulse solution is 
spectrally stable. In this case the zero A of £/o cancels the pole of £ v ,o at A. 


So, besides the position of the zeros of Syp and £ v o with respect to the imaginary axis, one 
should determine whether zero-pole cancelation occurs at the points in Afffiyo) with positive real 
part in order to decide upon spectral (in)stability questions. The next proposition focuses on this 
issue. 


Proposition 3.16. Let A 0 be a simple zero of Gyp. Suppose there exists do > 0 such that 

B(A o ,6 0 )n |J N(G S , o(-,y)) = 0- 

yeS 1 


Then, there exists <5i > 0 such that for every 0 < 6 < S\, provided e > 0 is sufficiently small (with 
bound dependent on 5), either (i.) B(A,,, 6) contains precisely one y-eigenvalue for every y € S 1 
or (ii.) B(A 0 ,6) contains no spectrum of f £ . The prior (i.) is the case if 



d v H 2 (uo,Vb{z))ip^,i(z)dz = 0, 


(3.10) 


where (pojx) = (tp 2o j(x), tp.i^x)) is a non-trivial bounded solution to (3.5) at A = A 0 , or 



*f'A«,2(z)d u G(u 0 ,v h (z),0)dz = 0, 


(3.11) 


where t/M 0 (x) = (x), i/%, 2 W) denotes a non-trivial bounded solution to the adjoint equation 

of ( | 3. 5[ ) at A = To. 


The orthogonality relations ( 3.10| ) and ( |3. 1 1 [ ) imply that there is no zero-pole cancelation. How¬ 
ever, the converse is not true as pointed out in Section |4.1| One can show that the integrals in 


the right hand sides of (3.10) and ( |3.11[ ) appear as one of multiple factors in the principle part 
of the Laurent expansion of £^o(-,y). Although it is possible to write down the singular part of 
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the Laurent series of £s,o(-,y) explicitly at a zero A e N(Sjfl) (of higher multiplicity), we de¬ 
cide to postpone this to Section [6.3.3 for the benefit of exposition, since the involved expressions 
are rather complex (except in the case m = 1 - see Proposition |4.6[ ). Eventually, these singular 
parts provide a tool to determine precisely whether zero-pole cancelation occurs or not. Therefore, 
Proposition |3.16| is weaker - but better digestible - than the statements in Section 6.3.3 


3.4 Small spectrum 


The derivative of the homoclinic (vh(x), q^{x)) in \2.1\ is a non-trivial bounded solution of (3.51 
for A = 0 due to translation invariance. So, by Proposition|3.7|it holds £/,o(0) = 0. Since, we have 


f 


d v H 2 (uo, v h (z))v' h (z)dz = 0, 


there occurs no zero-pole cancelation at A = 0 by Proposition |3.16| This has yet been predicted in 
Remark [3T4] by the fact £ e (0,1) = 0. Hence, a separate study of the spectrum close to 0 - the so- 
called small spectrum - is necessary in order to decide upon spectral (in)stability. The following 


Corollary of Theorems 3.10 and 3.11 makes this statement rigorous. 


Corollary 3.17 (Spectral stability is determined by the fine structure of the small spectrum). Sup¬ 
pose the following conditions are met: 

1. There exists 6q > 0 such that 

IJ M£*,o(-,y)) c {A € C : Re(T) < -<5 0 }; 

yeS 1 


2. For all zeros A 0 € A/f£/,o)\{0} of £/,o with Re(d 0 ) > 0, £ v ,o(-, y) has a pole of order 
oi'd,i, (£/,o) at A<> for every y e .S' 1 . 

Then, there exists c)] > 0, such that for all 0 < 5 < <5], provided £> 0 is sufficiently small (with 
bound dependent on 5), it holds 

cr(£ e ) nfleC: Re(T) > 0} c B( 0, 6). 

Moreover, there are for every y e .S' 1 at most ordoffiyo) y-eigenvalues inside /?((), 5). 


Remark 3.18 ( Determining the fine structure of the small spectrum ). In B71 the fine structure 
of the small spectrum is unraveled for nearly homoclinic wave trains. An expansion of the circle 
of critical y-eigenvalues around 0 is given in terms of their period. In our case only the (v, q)- 
components of the solution are close to the homoclinic (lyCr), qfix)) by assumption (S3) and 
Remark |231 Therefore, the results in ll47l are not directly applicable to our situation. Furthermore, 
in the Gierer-Meinhardt regime the structure of the small spectrum is known If53l Section 5]. 
However, the construction in ll53l depends essentially on the fact that the slow limit problem is 
linear and the model has two components. Hence, a direct extension to our (slowly nonlinear) 
situation is not immediately clear. This is subject of future research (Section |T2j). ■ 


Remark 3.19 ( Nonlinear stability). Since the spectrum cr(£ e ) always touches the origin due to 
translational invariance, nonlinear stability of is a delicate issue. On the one hand, one can 
‘factor out’ the translational invariance by looking at orbital (in)stability l35l . Then one measures 
the distance from the perturbation to the family of all translates of the solution rather than to the 
solution itself. On the other hand, one could look at diffusive nonlinear stability Il48ll50l . In that 
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case the perturbations are spatially localized, i.e. they are in certain Sobolev spaces, where the 
norms are adjusted by an appropriate weight p(x) = (1 - x 2 )' 1 for some d > 1/2. This allows for 
purely diffusive behavior of the perturbed solution. 


The presence of spectrum cr{£, £ ) in the positive half-plane jd e C : Re(d) > 0} implies orbital 
instability by 031 Theorem 4.3]. Therefore, the spectral instability criteria in Corollary |3. 12 


im¬ 


ply nonlinear orbital instability immediately. On the other hand, if one assumes, in addition to 
spectral stability, that 0 is a simple zero of £yo and that the spectrum <r(£., : ) touches the origin 
in a quadratic tangency, then we achieve nonlinear diffusive stability by If48l Theorem 1] or ll50l 
Theorem 1.1]. Therefore, also nonlinear stability is eventually determined by the fine structure of 
the small spectrum. ■ 


4 The cases n - 1 and m = 1 

Most specific models of the form <o> in literature have n = I or m = 1 so that the limit problems 
( |3.5| ), ( |3.6[ ) or (3.7 1 become 2-dimensional. In either case a more refined analysis of the fast and 
slow reduced Evans functions Gyp and £ v q is possible and allows for further insight into zero-pole 
cancelation and instability criteria. In this section we assume without loss of generality D\ = 1 in 


the case m = 1 and Di = 1 in the case n = 1 - see Remark 2.1 We start with an analysis of the 
fast reduced Evans function for n = 1. 


Proposition 4.1 (The case n = 1). Suppose n = I. All zeros of the fast reduced Evans func¬ 
tion £/.o: C,\ —» C are real and simple. There is precisely one positive zero do o/£/,o- The 
v-component of any non-trivial, bounded solution to (13.51) at A = do has no zeros. Moreover, 0 is 


a zero of&fp and (v’fx), cffx)) spans the space of bounded solutions to (3.5) at A = 0. 


Proof. We refer to Theorem 6.4 for the fact that every non-trivial, bounded solution to (3.5 l for 


d e Ca is a solution to the Sturm-Liouville eigenvalue problem. 


dG 


d xx v = (—(wo, i’h(x), 0) + A | v, v e H Z (R, ( 


(4.1) 


By |[32l Theorem 2.3.3] the set of eigenvalues d^ < ... < do of ( |4. 1 [ ) is finite, real and strictly 
decreasing. Moreover, the eigenfunction corresponding to d;, i = 0, ..., N is unique up to scalar 
multiples and has precisely i zeros. Hence, all zeros of £jo are real and simple by Proposition |3 .7 1 
Furthermore, the derivative v'f x) is a solution for d = 0 of (4.1). The function v’fx) has precisely 
one zero. So, we derive di = 0 and do > 0. □ 

We shift our focus to the case m - 1. Then, the slow reduced Evans function £ Vi o(d, y) is a 
quadratic polynomial in y, which gives rise to the following result. 

Proposition 4.2 (The case m =1). Suppose m — 1. The slow reduced Evans function £ v q : [Ca\ 
N(£>ffl) | X C —» C is given by 

fis.oCd, y) = y 2 - Tr(T(A)T s (2Lo, 0, A))y + 1. 

In particular, A e Ca\A((£/.o) is a zero of E> s p(-, y)for some y e .S' 1 if and only if 

t(A) := Tr(T(A)T s (2Lo, 0, A)) e [-2,2], (4.2) 


Proof. This follows readily by expanding £. v .o(d, y) as quadratic polynomial in y. 
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Example 4.3 ( Generalized Gierer-Meinhardt equation ). In 1531 the spectral stability of spatially 


periodic pulse patterns is studied, where (2.3) is the slowly linear generalized Gierer-Meinhardt 
equation ( 2.1 3| ) with f(u) = qu for some q > 0. Thus, system ( |3.7| ) is autonomous. Condition 
( |4.2| ) simplifies in that case to 

2 cosh<2£|) V7T2) + gWsi-hpr^FTl) £ ( _ 2 _ 2]> 

xfq + A 


where Q(A) is defined in (3.8). Although derived with a different method, this result agrees with 
Theorem 1.1.1]. ■ 


Recall that (2.4 1 is A'-rcvcrsible by Remark 2.2 Suppose 0 p£ (L e ) is contained in the space ker(7 


R ), whic h is th e case for the solutions considered in e xiste nce result Theorem 2.1 1| Combining 
this with (S3) -3 and the Hamiltonian nature of system (2.6) yie lds t hat u s is symmetric about Lo, 
i.e. it holds u s (Lq + x) = u s (Lq - x) for all iel. Hence, system (3.71 obeys a reversible symmetry. 
Thus, if <q{x, A) is a solution to (3.7), then so is x R s <p(2Lo -x, A), where R s is defined in Remark 
|2.2[ Hence, there exists non-trivial solutions u+(x, A) and n_( x, A) to 


drxU = (d u H i (u s (x), 0,0) + A) u. 


(4.3) 


which are symmetric and antisymmetric about Lo, respectively. This allows for further simplifica¬ 
tion of condition (4.21 in Proposition |4.2| 


Proposition 4.4 (The case m = 1 with reversible symmetry). Suppose m = 1 and u s is sym metri c 
about L 0 . Take A 6 Ca\A7(£/,o)- Let u+(x, A) and u~(x,A) be non-trivial solutions to (4.3), 


which are symmetric and antisymmetric about Lo, respectively. Condition (4.2) in Proposition 
simplifies to 


4.2 


Tr(T(T)7“j(2Lo, 0, A)) = 


— {u + (x, A)u-(x, T)](0) - @(A)u + (0, A)u-{ 0, A) 
dx 


AVIA) 

where r lV(A) denotes the Wronskian ofu+(-. A) and ;<_(•, A) and 

8H 2 


€[-2,2] (4.4) 


0(X) = 


f 


du 


■(n 0 , v h (x)) + d v H 2 (iio, v h (x))v in (x, A) 


dx. 


with vi n (x. A) the unique (meromorphic) solution to the inhomogeneous problem, 

dCj 

D 2 d xx v = ( d v G(u 0 , v h (x), 0) + A) v + —(n 0 , v h (x), 0), v e C". 

on 


(4.5) 


Proof. Express T s (2Lo, 0, A) in terms of the symmetric and antisymmetric solutions and invoke 
Proposition |4. 2 [ □ 

Example 4.5 ( Sinusoidal Gierer-Meinhardt model). In the companion paper ifTOl a sinusoidal 


Gierer-Meinhardt equation (2.13) with f(u) = sin(n) is studied. In this special case (4.3) is known 
as Lame’s equation. The symmetric and antisymmetric solutions u± to this equation are known 


in terms of Jacobi Theta and Zeta functions. Moreover, the homogeneous fast limit problem (3.5) 


is solvable in terms of hypergeometric functions. Applying the variation of constants formula 
gives the solution v,„(x, A) to problem ( |3.6[ ). Therefore, all terms in the right hand side of ( |4.4[ ) are 
explicitly known. This example shows that also in the slowly nonlinear situation the spectral sta¬ 
bility analysis can be very explicit. The sinusoidal Gierer-Meinhardt equation serves as a guiding 
example in ifTOll to study the spectral phenomena occurring when the periodic patterns approach 
the homoclinic limit - see also Section 17721 ■ 
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As mentioned in Section |3.3[ it is possible to obtain explicit expressions of the principle part of 
the Laurent series of <5 iSj o(-, y) at a zero d e N(Sf,o). Because of their complexity the expansions 


are treated separately in Section 6.3.3 However, in the case m = 1, the expressions simplify 


significantly. Therefore, it is worthwhile to devote a separate Proposition to this case. 


Proposition 4.6 (Singular part in the case m = 1). Suppose m = 1. Let d<> be a simple zero of&f$ 
and y € S *. The singular part of the Laurent expansion of£ s o(A, y) at A = A 0 is given by 


-yu(2Lp, dp) 
A - A<> 


r°° r°° „ 8G 

o v H 2 (uo, v h (x))v A Jx)dx Vy(x) —(uo,v h (x),0)dx, 

—CO —CO 


where u(x, A 0 ) is the solution to ( 4.3) at A = do with initial values u( 0) = 0, u'( 0) = 1 and vy is a 
normalized bounded solution (having L 2 -norm l) to 


Did xxV = (d v G(u Q , Vh(.r), 0) + do) v, v € 


(4.6) 


Proof. The statement is proven in a more general setting in Theorem 6.8 4. 


4.1 Robustness of zero-pole cancelation and an example for non-cancelation 

Proposition |4.6| shows that for m = 1 the slow reduced Evans function £ Vi o(-,y) has a removable 
singularity at a simple zero A,, of <3/ p if a nd only if one of the identities (3.101. (3.11) holds true or 
there exists a non-trivial solution to (4.31 at d = A 0 with boundary values u( 0) = 0 = li(2Lq). The 
set of d 0 € C for which (3.10) or ( |3. 1 I [ ) holds true will in general be discrete, since the involved 
expressions are analytic in d 0 . Moreover, lf60l Theorem 4.3.1-6] shows that this is also the case 
for the set of d<> e C for which the boundary value problem ( |4.3[ ), w(0) = 0 = u(2Lq) admits a non¬ 
trivial solution. Hence, zero-pole cancelation is a robust phenomenon in the absence of additional 
structure (such as the translational invariance at d = 0 mentioned in Remark[3.14[). 


Being robust, zero-pole cancelation can still fail in one-parameter families. Consider the Gierer- 
Meinhardt equation ( |2.13 ), where ai + 0 and f(u) = -pu with p > 0. We emphasize that in 
this case the slow reduced system (T6) is linear of center type. This differs from the ‘standard’ 
Gierer-Meinhardt setting considered in lfl2ll28lH6ll53ll59l . where f(u) = +pu and the slow re¬ 
duced system is linear of saddle type. 


Let «o > 0- Note that ( |2.13[ ) satisfies |(S1)H(S2)| and |(E1)| with Vh(x, no) > 0 for all x € R. 
Take u\ < 0 such that pu 2 = fT(uo) 2 + pu 2 . One easily checks that the four conditions in Remark 
2.10 holds true. Therefore, assumption |(E2)| is satisfied. Hence, Theorem |2. 1 1 1 implies that, for 
e > 0 sufficiently small, there exists a 2L e -periodic solution to (2.1 3J) satisfying assertions 
(S3)} 1,2,3. Moreover, it holds by Remark 12.121 


Lq(m) = ^ + sin_ 


mo 


V 


J(uq) z 


+ Ui 


In ll53l Lemma 3.3] it is shown that do = 1 /4(/L +1) 2 - 1 > 0 is the positive zero of the fast reduced 
Evans function &f,o- Note that both dyHiiuo, i’h(v)) and d u G(iio. vfx), 0) are strictly negative for 
all x e R. Moreover, the v-component of any non-trivial solution to (3.5 1 at d = do has no zeros 
by Proposition |4.1 1 Therefore, identities ( 3. 10[ ) and ( |3. 1 1 } are not satisfied. Now assume p > do- 
The solution to ( |4.3| ) at d = do with initial values u( 0) = 0, u'( 0) = 1 is given by, 

sin( yjp - Aqx) 


u(x, d 0 ) = 


VP - 4 0 
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Clearly, it holds u(2Lq, /l 0 ) = 0 if and only if 


/ kn \~ 

p = do + —- , 

\2Lo(n)) 


(4.7) 


for some k e Z>i. Since Loin) e (?r/2 ,tt) for every p > 0, equation (4.7) will have a solution 
p = Pk > 0 for every k e Z>i. We conclude with the aid of Propositions |3J6| and |4~6| that for every 
q > 0 there exists 5^ > 0 such that the disc B(Ao, S M ) contains, provided e > 0 is sufficiently small, 
either no spectrum of X e in the case n g [pk ■ k e Z>i) or a y-eigenvalue of for every y e .S' 1 in 
the case p - Pk for some & e Z>j. 


The transition of p through a point pk may seem like a blue sky catastrophe, which makes the 
pulse solution e ‘suddenly’ spectrally unstable. However, such a transition from cancelation to 
non-cancelation is caused by unstable spectrum moving through the point To- This can be seen 
by taking p sufficiently close to pk and comparing the values of the trace Tr(Y(T)7~ s (2Lo,0, T)) at 
A = p - [kn/ilLoin))] 2 and when T approaches To. The intermediate value Theorem and Proposi¬ 
tion [T6] imply that there must be a zero of £ lS ,oO, y) close to To for some y e S *. 


4.2 An instability criterion 

As a further application of our analysis, we show how to test for spectral instability using a parity 
argument. More precisely, one can gain insight in the values of f(T), defined in ( |4.2[ ), at a pole 
using Proposition |476] at T = 0 due to translational invariance and by taking the limit for T e R. to 
infinity. Subsequently, one determines with the mean value theorem, if the graph of t(A) crosses 
[-2,2], The use of such a parity argument is common in Evans function based stability analyses. 
See l45l Chapter 4.2] for an overview of the literature on this topic. 


Corollary 4.7. Suppose m = n — 1. Let To be as in Proposition 4.1 Denote by u(x, To) the 
solution to (4.3) at A = To with initial values u{ 0) = 0, //(()) = 1. Moreover, let v^ 0 be the 
normalized bounded solution (having L 2 -norm 1 ) to (4.6 1 at A = To- Finally, let \’h(x, u)for u € Uh 
be as in Remark [2777] Consider the quantities, 


1 1 := m(2Lq 


,*>) f 

—o 


oh 2 

If 


(uo,v h (x))v Ao (x)dx 


X O 

o 


dG 

du 


(uo,v h (x),0)vMx)dx, 


1 2 u'(0)(C?()u'(0) - 2H\(u s (0) 


,0,0)) f 

Jo 


,Zo (d u H x (u s (x),0,Qi) + 1 )[(u'(x)) 2 - (H^ixlOM 2 ] 


[(u' s (x)) 2 + (HiiusixlO^)) 2 ] 2 


dx 


+ 


where 


C 2 + QqHx (« s (0), 0,0K(0) - {Hi (m s (0), 0,0)) 2 
( m '(0)) 2 + (7/i(m s (0),0,0)) 2 

8H 2 


Go = 


-f 


du 


(u, v h (.r, u)) 


dx. 


U=UQ 


Then, the periodic pulse solution 0 p-e is spectrally unstable, for e > 0 sufficiently small, if: 
(i.) 1 1 <0 or (ii.) u s is symmetric about Lq and it holds 1 2 > -1. 


Proof. By Proposition |T2] the slow reduced Evans function G x p(A, y) has a zero for some y e .S’ 1 
if and only if fl4.2| ) holds true. Our approach is to gain information about the trace t(A) at T = 0, 
T = To > 0 and when Tel tends to infinity. An application of the intermediate value Theorem 
will eventually yield the result. We distinguish three cases: (1.) I\ < 0, (2.) I\ = 0 and (3.) 
Jo>0, J2> _ 1 and u s symmetric about Lq. 
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The case I\ < 0 

By combining Propositions ^. 2| and |4.6| the singular part of the trace t(A) at A = do £ Af(£/,o) equals 
I\. On the other hand, it is shown in Lemma |C2] that t(A) tends to infinity as A —> oo. Therefore, 
Proposition |4. 1 | and the intermediate value Theorem imply that there exists d* e (do, oo) such that 
t(d*) e [-2,2], if I\ <0 holds true. Hence, the periodic pulse solution 0 p e is, provided e > 0 is 
sufficiently small, spectrally unstable by Proposition [4T2] and Corollary |3.12|-1. 


The case I\ = 0 

If we have 1 1 = 0, then j) has no pole at do £ AA(£/,o) for each y € S 1 by Proposition |4.6| 
which implies is spectrally unstable, for s > 0 sufficiently small, by Corollary |3.12[-2. 


The case I\ > 0, i~ 2 > -1 and u s symmetric about Lq 

At d = 0 the derivative m'(x) is due to translational invariance a non-trivial solution to (14.3b, which 


is anti-symmetric about Lq. By Rofe-Beketov’s formula ||4] Chapter 1.9] a 1 inear independent 


solution to (4.3 1 is given by 


z(x) := w'(x) I 
JLo 


(d „//i(u s (y),0,0) + l)[«(v)) 2 ~ (//i(u s (y),0,0)) 2 ] 
[(«'Cv)) 2 + (^i(»sCv),0,0)) 2 ] 2 
Hi(u s (x), 0,0) 


dy 


(m'(x)) 2 + (//i(« s (x),0,0)) 2 


Note that z(x) is symmetric about Lq and that the Wronskian of z and u' s is constant with value 1. 
Furthermore, there exists C e R. such that the solution v,„(x, 0) to (|4.5[) is given by 


Vm(x, 0) = (x, u) 

ou 


+ Cv' h (x). 


U=UQ 


Substituting the obtained identities in the expression for the trace in Proposition |4,4| yields /(()) = 
21 2 - Again, the intermediate value Theorem implies that there exists a d* € (0, do) such that f(d*) £ 
[-2,2] if 1 1 >0 and 1 2 > -1 holds true. Hence, the periodic pulse solution d) pj - is, provided e > 0 
is sufficiently small, spectrally unstable by combining Proposition |4.2|and Corollary |3.12|-1. □ 


5 The Riccati transformation 


The Riccati transformation provides a tool for the decoupling of singularly perturbed linear sys¬ 
tems by block diagonalization. The decoupling takes place via a linear non-autonomous transfor¬ 
mation, which is determined by two matrix Riccati equations. Specifically, for bounded matrix 
valued A,y(x, e) the coupled system. 


( d x <p = e(A\i(x, e)ip + A i2 (x, e)ip) + „ 2 

W = A 21 (x, e)ip + A 2 2 (x, e)if/ ’ 

with 0 < e 1, decouples via a coordinate change (the Riccati transformation) to 

( dxX = e [Ah(x, e) + Aj 2 (x, e)U e (x)\x , , c „ 1+ „ 2 

\d x m - [A 22 (x,e) - eU £ {x)A l2 (x,e)]a) 


(5.1) 


(5.2) 


where U f (x) is a family of matrices satisfying a certain matrix Riccati equation as detailed be¬ 
low. Applying this decoupling to the linear stability problem ( |3.3| ) is a major step in reducing our 
stability problem to a ‘fast’ and a ‘slow’ stability problem yielding the factorization ( |1.2| ) of the 
Evans function. 
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Although the construction of the transformation is based on two results of Chang |(7] Theorem 1] 
and f8j Lemma 1], the assumptions on the coefficient matrices in |[8j are too restrictive. Therefore, 
we need a refinement of his statements. For this reason and the fact that the Riccati transformation 
lies at the core of our analytic factorization method, we present the full construction of the trans¬ 
formation in the next theorem. Moreover, we also prove that a periodic coefficient matrix implies 
periodicity of the Riccati transform, which appears to be a new result - see Remark [572] 

Theorem 5.1. Let n \, n 2 e Z>o, to e R.>o andMj e Q>(R x (0, eo)> Mat„ jXnj .(C)), i, j = 1,2. Assume 
that 


d x ip - A 2 fix, e)f, f £ C" 2 , 


(5.3) 


admits an exponential dichotomy on R with constants Kq,pq > 0, independent of e. Then, for 
e > 0 sufficiently small, there exists continuously differentiable matrix functions Ufix') and S fix) 
satisfying the matrix Riccati equations, 


d x U = A 22 U - eUA n - eUA n U + A 2 u 
d x S = e(Au+A n U)S-S(A 22 -eUAn)-A n , 


l /22X«i (C), 


U £ Mat,, 

S £ Mat„ lX ;; 2 (C), 


(5.4) 


with the following properties: 

1. If- and S f - have e-independent bounds: 

\\U e \\ < 8K 0 p^\\A 2l \\, ||5 e || < 8 C//q 1 HA 12 II 
where C > 0 is a constant depending on Kq only. 

2. The coordinate transform, 

' X 


(5.5) 


= Hfx) 


LO 


, H e {x):= 


I -eSfx) 

Ufx) I - eU f (x)S f-(x) 


(5.6) 


diagonalizes system (5.1) into (5.2). 

3. Let 9 £ (0,1). The unique bounded solution Q, to the inhomogeneous matrix problem. 

d x D.-A 22 (x,e)Ll+A 2 i(x,e), Q £ Mat„ 2Xni (R, C), (5.7) 

satisfies for x £ R 

\\U £ (x)-n £ (x)\\<Ce l ~ e , 

where C > 0 is a constant depending on Kq,pq, ||An||, ||A 12 II and HA 21 1| only. 

4. Let a > 0. For each x £ R, we have the approximation, 


11140011 < — 

Fo 


sup (e||A 1 2lll|l/ e 0v)ll 2 + l|A 2 tCy, e)ll) + Ce~ a ^ 2 

ye[x-a,x+a\ 


(5.8) 


where C > 0 is a constant depending on Kq,pq, 11A 12 11 and HA 21 1| only. 


5. If the matrices Aij(-,e) are L-periodic for 1 < i,j < 2, then the coordinate transform H ( is 
also L-periodic. 
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Proof. First, we set up an integral equation for U e and prove global existence via a contraction 
argument. Since U e triangulizes the system, an integral equation for S e can be derived from the 
variation of constants formula. The first four properties of U e and S f follow readily from the 
integral equations they satisfy. Finally, periodicity of the transform is proven by exponential sep¬ 
aration. 

Existence of U f -(x). 

Since An is uniformly bounded in e on R, system. 


d x (f = eA 11 (x, e)<p, ip e C" 1 , 


(5.9) 


has bounded growth on R with constants K\,ep\ >0, where K\ = I and p\ is independent of e. 
Denote by T\ e {x,y) and T 2e {x,y) the evolution operators of systems (5.9) and ( |5.3| ), respectively. 
Take p = SK^p 0 1 IIA 21 1|- The ball B((),p) c C/,(R, Mat„ 2XH] (C)) is a metric space endowed with the 
supremum norm. We want to show that the map : B( 0, p) —> B( 0, p) given by 


Oft 


■U)(x)= f T£ e (x,y)[-eU(y)A 12 (y,e)U(y)+ A 21 (y,e)\T li£ (y,x)dy 

— CO 

r* 00 

- I Tl(x,y) [~eU(y)A l2 (y, e)U(y) + A 21 (y, e)] T u (y, x)dy, 

J X 


is a well-defined contraction. If e > 0 is sufficiently small, it holds for all U e B( 0, p) 

2K 0 


mu\\ < [6P 2 ||A 12 || + ||A 21 ||] < p. 


Po - epi 

Therefore, J\ f is well-defined. Similarly, provided e > 0 is sufficiently small, we estimate for 
UuU 2 eB(0,p ) 

\\JK e U x -JK e UA < 4e *° Pl|Al2ll ||£/i - U 2 1| < ||I/i - U 2 \\. 

Bi ~ epi 

Hence, J\ e is a contraction mapping. By the Banach fixed point Theorem the integral equation 
cAJJ = U has a unique solution U f -(x) in B(0,p). It is readily seen by differentiating this integral 
equation that U e satisfies the matrix Riccati equation ( |5.4[ ). Moreover, we derive the bound on U e 
in 


. Since (5.3) has an exponential dichotomy on 1, (5.7) admits a unique bounded solution 
by Proposition B.12 By Proposition [ bT 2J it holds for \x -y\ < e~ e 


\\T u {x,y)-I\\<C'e l - e , 


(5.10) 


where C' > 0 depends on ||An|| only. Using (5.10) and the bound on U e in (5.5 1 , we estimate for 

rel 


\\u e (x)-n e {x)\\ = 


X x r* 00 

T s l e (x,y)A 2 i(y, e)dy + I T“Jx,y)A 2l (y, e)dy 

OO X 

4K 0 p- 1 (ep 2 \\A l2 \\ + 6 1 - e C'||A 21 || + e -^|A 21 ||). 


This shows the third property. The fourth property follows by splitting the interval of integration 
of the two integrals in the right hand side of the identity U e {x) = (jA e U e )(x). We obtain four inte¬ 
grals over (- 00 , x - a), (x - a, x), (x, x + a) and (x + a, 00 ), respectively. We estimate each integral 
separately and obtain approximation (5.8 1 with C = p 2 ||Ai 2 || + ||A 2 i||. 
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Existence of S £ (x). 

Since An, An and U £ are uniformly bounded in e on R, system, 

d x x = e[A n (x,e)+An(x,e)U e (x)\x, feC" 1 , (5.11) 


has bounded growth on R with constants K 2 , ep 2 > 0, where K 2 — 1 and p 2 is independent of e. 
On the other hand, equation, 


d x o = (A 22 (x, e) - eU £ (x)A\ 2 (x, e))a>, co e C" 2 , 


(5.12) 


can be seen as a perturbation of (5.3). By roughness (Proposition |B.9| ) it therefore possesses an 
exponential dichotomy on R with constants K^,p 4 > 0. Here, K 4 depends on Kq only and we 
choose IJ 4 = po/2. Denote by T 2 e (x,y) and T^ e (x,y) the evolution operators of systems ( 5.1 1 [ ) and 
(|5 .1 2[), respectively. We define .S' f (x) via the variation of constants formula, 


X X S' OO 

T^ e {x,y)A n (y,e)TlJy,x)dy + I T 3>e (x, y)A , 2 (y, e)T s Ae (y, x)dy, 

OO %J x 


We readily derive the bound on S £ in (5.5 1 . This proves the first property. It is easily verified by 
differentiation that S £ satisfies the matrix Riccati equation (|5.4[). Finally, using S £ and U e satisfy 


equations (5.4), it is a straightforward calculation to see the change of variables (5.6) transforms 
system (5.1) into (|5.2[). This proves the second property. 


Exponential separation and periodicity 

Only the fifth property remains to be proven. Our plan is to show that system <HD is exponentially 
separated in the sense of ||39ll . Subsequently, we make use of the fact that exponential separation 
preserves periodicity. Therefore, denote by P £ (x), x e R the projections corresponding to the ex¬ 
ponential dichotomy of system ( 5.1 2[ ) on R, established in the latter paragraph. We define the 
following projections, 


P i, e = 


0 0 
0 Pe( 0) 


P 2, e = 


I 0 
0 0 


P\e = 


0 


0 


0 l-Pel 0) 


Denote by f c C" 1+ ” 2 the range of the projection P l£ for i - 1,2,3. Let nn be the rank of P e (0). 
Using both the bounded growth of ( [5.11[ ) and the exponential dichotomy of (5.12), we conclude 
system that ( |5.2| ) is, for e > 0 sufficiently small, (nn, n\, in - m 2 )-exponentially separated with 
respect to the decomposition V\ t£ © VS.e © Vye As a result, system (5.1) is also (m 2 , n\, « 2 _ m i)- 
exponentially separated with respect to the decomposition W i, e © Wo.e © VTi f , where W Lr is the 
range of the projection Qi £ := H £ (0)Pj^ £ H~ l (0) for i = 1,2,3. 


Now, suppose Ajjf, e) are L-periodic for 1 < i, j < 2. Let X t (x) be the principle fundamental 


matrix of system (5.1) with X e (0) = 7. Invoking (5] Corollary 4] gives that X £ (-)Q 2 f X f (-) 1 is 


L-periodic. Denote by T £ (x,y) the evolution of the diagonal system (5.2). We calculate for x e 


X £ (x)Q 2 MxY X = H £ (x)T £ (x,0)P 2 , £ T £ (0,x)H £ (x)-' = H £ (x)P 2 , £ H £ (xy l 

1 - eS £ (x)U e (x) e,S £ (x) 

U e (x) + eU £ (x)S £ (x)U £ (x) eU £ (x)S £ (x) 


Hence, S £ , U £ S £ ,S £ U £ and U £ + eU £ S £ U £ are L-periodic. So, U £ S £ U £ is also L-periodic. Combin¬ 
ing this with the L-periodicity of U £ + eU £ S £ U £ , we conclude that U £ is L-periodic. This implies 
that H £ is L-periodic, which concludes the proof of the fifth statement. □ 
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Remark 5.2 ( Periodicity of the transform). The periodicity of the transform in Theorem |5.1[ -5 
is a new discovery to the authors’ knowledge. It is natural to ask whether there always exists a 
periodic choice for a coordinate change, which transforms a periodic system into diagonal form. 
However, it is shown in ll38l Chapter 5] that this is not the case. It seems that the periodicity of 
the coordinate change H e is due to the special (slow-fast) structure of system ( |5. 1 [ ) used in the 
(non-generic) proof of Theorem |5.1|-5. ■ 


Remark 5.3 ( Exponential trichotomies). The concept of exponential dichotomies can be gener¬ 
alized to exponential trichotomies iBTI . In addition to a stable and unstable direction, a center 
direction is considered. In fact, the {m 2 ,n\,n 2 - «? 2 )-exponential separation of systems ( |5.1[ ) and 
(5.2 1 obtained in Theorem 5.1 shows that these systems admit an exponential trichotomy on R. 
Here, the stable and unstable directions are related to the fast subsystem and the center direction 
is related to the slow subsystem. ■ 


Remark 5.4 ( The Riccati transform in a broader setting). Formally, the Riccati transform can be 
employed to diagonalize general linear equations as pointed out in El Remark 4.7]. However, 
the Riccati solutions can become singular in finite time. We use both the slow-fast structure of 
( |5.1| ) and the exponential dichotomy of (5.3 1 to achieve global boundedness of the transformation 
functions U e and Sm 


6 Proofs of the main results via the analytic factorization method 


6.1 Approach 

As mentioned in the introduction we employ an analytic factorization method to prove our main 
Theorems 3.10| 3.11 and 3.16 We start by showing that the spectrum of X £ is contained in an 
e-independent sector. This provides an important a priori bound on the magnitude of the spectrum. 


The construction of the e-independent sector is carried out in Section 6.2 


Subsequently, we start constructing the reduced Evans function £o: C,,\ xC -» C, whose zero 
set approximates the zero set of the Evans function £ £ . As stated in section |3.2| £(> is defined 
in terms of the three singular limit problems (|3.5|), (3.6 1 and (|3.7[) of the linear stability problem 


(3.31. We study these three singular limit problems first separately in Section [63] before drawing 
the connection between £q and the Evans function £ e . 


The reduced Evans function £o is a product (3.4) of an analytic fast reduced Evans function £y o 
and a meromorphic slow reduced Evans function fi v o. This factorization of £o corresponds to 
the fact that system 


can be diagonalized using the Riccati transform. Indeed system (3.3) is 
clearly of the form ( |5.1[ ). However, the application of the Riccati transformation is only legitimate 
in the case system, 

■n2« 


d x f = ft22A x > e 


( 6 . 1 ) 


has an exponential dichotomy on R. Using roughness techniques we can show that the latter is 
the case, whenever A is away from the discrete set of eigenvalues of the homogeneous fast limit 


problem (3.5), which correspond to the zeros of the fast reduced Evans function £y_o- This result 
is the content of Section 6.4 Consequently, using the periodicity of the coefficient matrices of 
system (3.3), it is possible to factorize the Evans function £ £ into two factor £/ /; and £ vg . This is 
performed in Section 6.5 The two factors £/ j£ and S S E can be linked to S/,o and G s fl, respectively. 


The latter is the content of Section 6.6 Finally, an application of Rouche’s Theorem yields the 


proofs of Theorems 3.10 3.11 and 3.16 This can be found in Section 6.7 
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6.2 A priori bounds on the spectrum 

In this section we show that the spectrum of X £ is contained in an e-independent sector. 


Proposition 6.1. The operator £ E : C^(R, R"' + '') c C ub ( R,R m+n ) -» C ub ( R,R m+ ") given by 
is sectorial and densely defined. For e > 0 sufficiently small, there exists constants a> e R>o and 
m € (7r/2, n), both independent of s, such that the sector X := {A € C : A ± a>, |arg(d - oj)\ < 
vj) U {m} is contained in the resolvent set p(-C f: ). 


Proof. It follows by K33i Corollary 3.1.9.ii] that Jf E is densely defined. An application of 
Theorem 1.3.2] shows that X £ is sectorial as sum of a sectorial and a bounded operator. At first 
sight this seems to be sufficient. However, the matrix B f: defined in (3.2) will in general have a 
norm of order Oie '). Therefore, we need to follow a different path to prove that the sector corre¬ 
sponding to X e may be chosen independent of s. 


Our approach is to decompose X £ in more elementary building blocks to control the s 1 -terms 
in B e . First, we show that the operator Xi, £ : C 2 6 (R, R"') c C ub {R, R m ) —> C ub (R, R m ) given by 

£i, e O) = (D\d n + E~ l d u H 2 (up,s(-),Vp, E (-)))u, 


is sectorial with an e-independent sector. Subsequently, we prove this for £ e : C\ (R, R m+ ”) c 
C ub ( R, M m+M ) C ub ( R, R m+ ") given by 

/ u \ ( Xi, £ (m) + s- l d v H 2 (up, E (-), v p , £ (-))v \ 

\ v / ^ \ £ 2 D 2 dxxV / 1 

Finally, X e can be seen as a perturbation of L e by a bounded operator with Oi I )-norm. 


The operator Xi, £ 

Our goal is to show that the spectrum of Xi, £ is contained in an e-independent sector. By G3l 
Proposition 2.1] it is sufficient to show that the associated eigenvalue problem. 


f yfDidxU 
| sfDldxP 


^ | d “ H -^V-d x )Xy.dF)) j u » u,peC m , 


( 6 . 2 ) 


has no non-trivial bounded solutions for A in some e-independent sector. Denote by T\ f: (x, y, A) 


the evolution of system (6.2). Let T [ (x,y,A) be the evolution of 


fTfdfU = Vd/; 
y/DidxP = VAu 


U, p € 


(6.3) 


By solving system ( |6.3[ ) explicitly, one observes that, whenever A e C\R<o, system (6.3) has 
both bounded growth with constants K,px,+ > 0 and an exponential dichotomy with constants 
K,px- > 0 - Here, we have px.+ = ||D 1 _1 ^||JRe( Vd)|, p,\- = ||Di|p^ 2 |Re( Vfi)| and K - 1 . Since 
we have lim £ in L r = L 


-o 


such that for all A e S i := [p e 


> 0, as stated in |(S3) 


1, there exists a constant C\ > 0, independent of s, 
: |Re(-^u)| > Ci) it holds h\ \ = p~ x l sinh _1 (4) < L e . 


Let A € Si. Note that H 2 vanishes at v = 0 by |(S1)[ Combining this with (S3)| -3 yields K\,p\ >0 
such that 


\\duH 2 (up, s (x),Vp, s (x))\\ < K { e £ ^ lW 
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for x 6 [~L e ,L E ]. Let w, z € R such that 0 < z - w < 2h,\ < 2 L s . Taking into account the 
2L e -periodicity of n P;£ and v p e , we have 



ll^w-^2(Wp,e(x), Vp, £ (x))|| 
-—- ax < 

Vi4e 


2K\ 

p i VuT 


By Lemma |BT2| we conclude for \w - z| < 2 / 1,1 

2 k ( 2,K 

\\T i(2, w, A) -T 1 >£ (z, w,4)ll < -exp 2 y u, l , + /7, l +- 

Pi VUI \ Pi VUI 

Since n^+h* is independent of A, we derive that there exists a constant Co > 0 such that, whenever 
A € S 1 satisfies |4| > C 2 , then it holds for all |w - z| < 2hx 


\\Ti(z,w,X)-T u {z,w,X)\\ < 1. 


(6.4) 


Now let Si be an e-independent sector disjoint from B( 0, C2) U [C \ S 1 ] - see Figure [3a] For all 
f e Ip there are no non-trivial bounded solutions to ( |6.2[ ) by combining ( |6.4[ ) with Proposition 
B.10I So, by Il23l Proposition 2.1] the resolvent set p(JA\. e ) contains the e-independent sector Si. 


The operator X £ 

Consider the elliptic operator X 2 : C^(R, R") c C ub { R, R”) —» C ub { R, R' ! ) given by = 

D 2 d M v. Clearly, we have p{£, 2 ) = C\R<o 3 Si. For A 6 Si the operator on C M /,(R,R m+ ”) 
defined by 


a \ ( (£ 1 , 8 -X) x (u-s l d v H 2 (up tE (-),v P ' E (-))(£?£ 2 ~ 4) *(v)) 

v ) " [ { e 2 £ 2 ~ drVv) 

is an inverse of £ s - A. Therefore, the sector Si is contained in the resolvent set p{£ E ). 

Conclusion 

Define 


BbA*) ~ 
Let £ b y. C„i(R,R m+ ") -» 


d u H 1 (n p>£ (x), v P;£ (x), s) d v H\ (u p , E (x), v p Jx), s) 
d u G(u p Jx), v P;£ (x), e) d v G(u p Jx), v p Jx), e) 


CubQ I 


l m+n ) be the multiplication operator [ £b, s T\(x) = &i,. E (x)<f- By 
Remark [23] the norm of L b ,s satisfies ||Xz», e || = 0(1). 


Invoking ll25l Theorem 1.3.2] and its proof yields the conclusion: the sum £ E = £ E + X/>, £ with 
domain C“, ( R, R m+n ) is densely defined and sectorial with an £-independent sector S c p{£ e ). 
Here, we have used ||X& )£ || = (3(1) and that Si is independent of s. □ 

Proposition |6.1| shows that the spectral stability of the periodic pulse solution 0 p e is fully deter¬ 
mined by the y-eigenvalues of £. E in one of the following family of regions in C - see Figure 

m 


Notation 6.2. Let a) e R>o and m e (tt/2, n) be as in Proposition [6T] For A e R<o we denote 
^A,o := {4 e C : Re(4) e (A, oS), |arg(/t - a>) \<n- m}. 
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(a) Construction of the sector E; in the proof 
of Proposition |6.1 1 



(b) The region E A ,o contains the spectrum 
o-(£.s) inC A . 


Figure 3: a priori bounds on the critical spectrum. 


6.3 The three singular limit problems 

As a prerequisite to properly define and analyze the reduced Evans function £o and its factors G/p, 


G s fl, we treat the homogeneous fast limit problem (3.5 l, the inhomogeneous fast limit problem 


(3.6 1 and the slow limit problem (|3.7|), consecutively. 


6.3.1 Homogeneous fast limit problem 


The homogeneous fast limit problem (3.51 can be seen as the singular limit of our linear stability 


problem (3.3), where the /(-component is taken constantly 0. To track solutions to (3.5) properly, 


we establish an exponential dichotomies of system (3.51 on both half lines. A sufficient condition 


for this is that the coefficient matrix of (3.5) is asymptotically hyperbolic. With the tools that 


have become available when (3.5) has exponential dichotomies on both half-lines, we are able to 
construct an analytic function Gf$ that detects the values of A e C for which (3.51 has bounded 


solutions. This will be the fast reduced Evans function. The above is the content of the following 
lemma and theorem. 


Lemma 6.3. There exists Ao > 0 such that for A e (—Ao, 0) the spectrum of the matrix, 

0 D~ ] 

d v G(u, 0,0) + A 6 


A(u,A) := 


is uniformly bounded away from the imaginary axis on Ku X Ca, where Ku is as in Remark 2.3 
and Ca as in Notation 


Proof. For a matrix A e Mat„ xn (C) denote by T(A) = {v*Av : v e C", ||v|| = 1} its held of 
values. Since d v G(u, 0,0) has positive definite real part by |(S2)[ the held of values K(d v G(u, 0,0)) 
is for every u 6 'Ku contained in the positive half plane by l26l Property 1.2.5a]. In fact, by 
compactness of ( Ku there exists Ao > 0 such that we have T(<fG(u, 0,0)) c C_a 0 f° r every 
u e 'Ku . Let A e (-Ao,0). For u e Ku and A e Ca we establish using |[26l Property 1.2.3] and 
lf26l Corollary 1.7.7] 

o-((d v G(u, 0,0) + A)D~ 2 l ) c (K(d v G(u,0,0)) + A)T{D 2 X ) c [peC: Ref//) > d^f A,> + A)}, 
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where d max is the largest diagonal value of D 2 . The eigenvalues of A(u, A) are given by the square 
roots of the eigenvalues of (d v G(u, 0,0) + A)D~ ] . Therefore, we obtain for u e 'Ku and A e C\ 
that any eigenvalue p £ <r(A(u, A)) satisfies |Re(ju)| > cos(7r/4) V(Ao + A)/d max , which concludes 
the proof. □ 


System Q has exponential 


Theorem 6.4. Let A € (-Ao,0) with Ao > 0 as in Lemma 6.3 
dichotomies on [0, oo) and (-oo,0 ] for every A £ C a with constants K r (A), p,(A) > 0, chosen 
continuously dependent on A. Moreover, one can choose the corresponding rank n projections 
r P r± (±x, A), x £ [0, oo) such that P r ±(x, •) is analytic on Ca- 


Let B l ‘: s : Ca —* Mat2 nX n(C) be analytic bases ofP r+ ( 0, A)\C 2n \ = B"(d)[C"] and ker(!P r _(0, A)) = 
B*(A)[C n ] for A e Ca- The analytic function Gfp : Ca —> C given by Gfp(A) - dct(B"(A). BfA)) 
has the following properties: 


1. £/;o(d) = 0 if and only if (3.5) has a non-trivial (exponentially localized) solution in 
C b ( R, C 2n ). 


2. Sffi(A) 4 0 if and only if ( |3.5|) has an exponential dichotomy on R. 


3. 

4. 


The zero set of N(Sfp) c Ca is discrete and independent of the choice of bases B l ‘ 


A zero Aof&fp is simple if and only if (3.5) admits a unique non-trivial solution in Cfo(R, C 2 ") 
up to scalar multiples. 


Proof. Let A e Ca- The matrix jdf () (A) := A(uo,A) is hyperbolic by Lemma [673 
: coefficient matrb 

\\^ 22 ,o(x, A ) - JK 


Recall that 


^ 22 , 0 (x,A) is the coefficient matrix of (3.51. By Remark |2.4| there exists K\,p\ > 0 such that 
r 220 (A)\\ < K\e~ Ml ^ for x £ 


. Following lf40l Lemma 3.4] and li46l Theorem 1] 
system p.5[) has exponential dichotomies on both [0, oo) and (-oo, 0] with the desired properties. 


By Il40l Proposition 2.1] we have £/,o(/l) 4 0 if and only if (3.5 I has an exponential dichotomy 
on R. On the other hand, every solution ip £ C&(R, C 2 ") to (3.51 must satisfy ip(6) £ B : fA)\C n \ n 


B;X/t)[C"]. This settles the first two properties. To prove the third property, we consider the oper- 


l n ) C*(R, R") given by 
£ r (v) = (D 2 d xx - d v G(u 0 ,v h (-),0))v. 


ator L r : C 2 (R, R M ) c C b (l 


By li25l Theorem 1.3.2] and 1331 Corollary 3.1.9.i] £ r is sectorial. The eigenvalues of the sectorial 
operator f r correspond to the zeros of the analytic map £y o- We deduce that fiy o has a finite num¬ 
ber of zeros in Ca- The third property has been proven. The fourth property follows immediately 
by EB Theorem 3.19]. □ 

Although the following corollary is an easy consequence of the boundedness of Xa,o> it will be of 
great importance to obtain d-uniform estimates later. 


Corollary 6.5. Let A £ (-Ao,0) with Ao > 0 as in Lemma 6.3 System (3.5) has exponential 
dichotomies on [0, oo) and (-oo, 0] for every A £ Ea.o with constants K r , p r > 0, independent of A. 


Another important consequence of the exponential dichotomies established in Theorem 6.4 is that 


the differential operator associated with (3.5 l is Fredholm. 


For each A £ Ca the bounded 
) —> C b ( R, C 2 ") given by £a(4 >) = ( d x - Jl 22 .o(-, Aj)p is Fredholm of index 


Corollary 6.6. Let A £ (-Ao, 0) with Ao > 0 as in Lemma 6.3 
operator Fa ■ C ; '(R, 

0. Moreover, Fa A invertible if and only if A £ C ; \ \ A'ffi/.o)- 
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Proof. This is the content of BTH Lemma 4.2], 


□ 


Finally, Theorem |6.4| provides us with the fast reduced Evans function and thereby proves Propo¬ 
sition [3/7] 


Definition 6.7. Let A e (-Ao, 0) with Ao > 0 as in Lemma 6.3 The map : C a —» C given by 
I'jjq(A) = det(/?"(/!), B s r (A)) obtained in Theorem 6.4 is called the fast reduced Evans function. 


6.3.2 Inhomogeneous fast limit problem 


The inhomogeneous fast limit problem (3.6 1 can be seen as a family of layer problems associated 
with the singular limit of the linear stability problem ( |3.3[ ). The matrix solution X in (x, A) to the 
inhomogeneous problem ( |3.6[ ) forms one of the key ingredients of the slow reduced Evans function 
S s fl- We will prove that A,•„(■*, •) is meromorphic for each x e R. Singularities of X- in (x, •) can only 
occur when A is one of the zeros of &f,o, which are obtained in Theorem 6.4 


Theorem 6.8. Let A € (-Ao,0) with Ao > 0 as in Lemma 6.3 There exists a unique continuous 
solution Xin : R X C\\ Af(£/,o) —> Mat2 nX 2m(C) to equation (3.6) with the following properties: 


1. Xj n (x, ■) is meromorphic on C\ for all x & R. 

2. Xi„(x, ■) is analytic on Ca\W(£/,o) for all rei 

3. Xi n (-,A ) is exponentially localized for each A e C\ \ AT£/,o), i.e. there exists constants 
Ki n (A),pi n (A) > 0, chosen continuously dependent on A, such that for all x £ R 

\\Xin(x,A)\\ < K in {A)e-^ m . 


4. Let do e C,\ he a simple zero of £/.o. Denote by <p,i (x) and fx, (x) non-trivial bounded 
(exponentially localized) solutions to ( |3.5| ) and its adjoint equation, 

d x f =(x,A)*, (A € Mati x „(C), (6.5) 


respectively, such that 



q j 4 >x,(z)dz = 1. 


There exists a neighborhood /? (, c Ca of A 0 and a mapping X,i j n : ^ Ma.t2^x2m(^X 

such that on R X B ^ we expand 

X in (x, A) = f fxSz)AA2\,o(z)dz + X A ^ in (x, T). 

A A o J-oo 


Here, X,i.j n (x, ■) is analytic on ffi for every x 6 R. Moreover, X A<> j n (-, T) is exponentially 
localized for every A 6 B\_. 


Proof. For A e Ca the operator Ex defined in Corollary |6.6| is Fredholm of index 0 and Ex is 
invertible if and only if A e C ; \\ATfi/.o). Therefore, E { 1 is analytic on C ; \\A , (£/.o) and meromor¬ 
phic on Ca by Il34l Theorem 1.3.1]. This settles the first two properties. 


First, (3.5 I has by Theorem 6.4 an exponential dichotomy on R for A e C ; \ \ ATfiyo) with con¬ 
stants chosen continuously dependent on A. Second, recall that G vanishes at v = 0 by |(Sl)| So, 
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by Remark [2~4] ^ 21,0 is exponentially localized. Combining these items with Proposition B.12 


establishes the third property. By (3] Theorem 3.19] A 0 is a simple eigenvalue of the operator 
pencil A £.a if and only if A 0 is a simple zero of Gyp. Hence, the fourth property follows by an 
application of Keldysh formula for XT 1 given in Il34l Theorem 1.6.5]. □ 


Remark 6.9 (Connection with Fredholm alternative). Let / e C/,(R, C 2n ). The Fredholm alterna¬ 
tive in KOl Lemma 4.2] states that the inhomogeneous equation, 

d x ip = &22,0( x ’ A )<P + /(*)> <P e C 2 ", 

has a bounded solution if and only if the solvability condition, 



\p{x)f{x)dx = 0, 


is satisfied for all bounded solutions ip to the adjoint equation ( |6.5[ ). This corresponds directly with 
the fact that X m {x, •) has a removable singularity at a simple zero To of Gyp if and only if we have 


by 


Theorem 6.8 4. 



^ 00 ^ 21,0 (z)dz = 0, 


It is possible to obtain expressions for the singular part of the Laurent series of X ln at a zero of Gf t 0 
of higher multiplicity by looking at a canonical system of generalized eigenvectors. However, for 
simplicity of exposition we treat only one special case in the Proposition below, which requires 
significantly less notation. The interested reader is referred to If34l Chapter 1] for the general 
set-up. 


Proposition 6.10. Let A £ (-Ao,0) with Ao > 0 as in Lemma 6.3 Suppose d v G{uo,v,0) is 
symmetric for all v £ V. Let A<> £ C\ be a zero of Gyp of multiplicity k Ar and let ipxpfx), 
..., (pk lt> p r fx) form a basis of the space of bounded (exponentially localized) solutions to (3.5) 
satisfying the orthogonality relations, 



Q 0 ) TjpSz)dz. = 6jj, 1 < i,j < k;y. 


There exists a neighborhood B y c Ca of A,, and a mapping X,yj n ■ R X B y —> Mat2„ X 2»4C), such 
that for each A £ B,y we expand 


Xin{tC, A) — X Ar J f! (.V, T) + 


A-A ( 


Go „o 

Yj VjaS x ) 

7=1 


TjMiz)* 


0 -I 
I 0 


^ 21,0 (z)dz. 


Here, X A<i j n (x, ■) is analytic on B A for every x £ R. Moreover, X,yj n (-, A) is exponentially local¬ 
ized. 


Proof. Let £.a be as defined in Corollary |6.6[ Since d v G(uo, Vh(x), 0) is symmetric for all x £ 


tp(x. A) is a (bounded) solution to (3.5 1 if and only if <p(x, A) T ( J ) is a (bounded) solution to (6.5). 

L, it follows <p(x, T) t = ip(x. A)* 


Moreover, because all solutions < p(x,A) of (3.5 1 are real for A £ 

for each x £ R and A £ C. On the other hand, we have by Ii40l Lemma 4.2] that / £ ran(X.i) if and 
only if 



P(z)f{z)dz = 0, 
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for all bounded solutions fix) to ( |6.5[ ). So, for a bounded solution ifi(x,A) = (v(x, A), q(x, A)) to 

f 


(3.5 i we have the implication, 
dJLx 


dA 


-<fi(x. A) = 


0 

v(x, d) 


€ ran(£,0 


l|v(z,d)|| dz = 0 


by pairing ■^£, / np(x, A) with the solution <p(x,A)*(j ( j j to the adjoint equation (6.5). Hence, 

^■tp(x,A) is contained in ran(£ 2 ) if and only if <p(x,A ) is trivial. So, the generalized eigenspace 
of Xx equals the eigenspace ker(X.i). With this knowledge the result follows from an application 
of Keldysh formula for Xj 1 given in If34l Theorem 1.6.5]. □ 

6.3.3 Slow limit problem 


Consider the slow reduced Evans function £< n - see Definition 3.8 In the next result we establish 


the properties mentioned in Proposition 3.9 In addition, we obtain explicit expressions of the 


singular part of the Laurent expansion of £^o close to a pole using Theorem 6.8 


Proposition 6.11. Let A e (-Ao,0) with Ao > 0 as in Lemma 6.3 The slow reduced Evans 
function S s q is well-defined and has the following properties: 

1. fi v o is analytic in both A and y on C,\ \ NiG/fi) X C and G Si o(-,y) A meromorphic on C\for 
each yeC. 

2. For every y e .S' 1 the map £ v o(-, y) is non-trivial. 

3. Suppose d<> is a simple zero of Gyp. Let <px 0 = (jp; p,i, </X, 2 )> = (<A/U,i> Ba. and 


X,i„,in as in Theorem 6.8-4. Define for A e B [ 

-r 

kJ —o 


<p ■= 

f- 


d v H 2 (u 0 ,v h (z))<PA <> Az)dz e C , 


^^a(z)d u G(uo, Vh(z), 0)dz e Mati xm (C), 


-f 

X oo 

[d u H 2 (u 0 ,v h (z)) + d v H 2 (u 0 ,v h (z)) r V , A»,in(z, A)]dz € Mat mXm (C), 

oo 


( 6 . 6 ) 


where denotes the upper-left (n x m)-block of the analytic (2 n x 2 m)-matrix Xx^in- 

Moreover, let (ufx. A), pfx, A)), i = 1,..., 2m be a fundamental set of solutions to the slow 
limit system (3.7). Finally, let C(d, y) be the cofactor matrix of 


dd a (A, y) : = 


I 0 

6a{A) I 


T s (2Lq, 0, d) - yl e Mat 2mX 2 m(C). 


(6.7) 


For all y € S , the singular part of the Laurent series ofG s p(-, y) at d 0 is given by 


1 2 

— X X)) (/ 

(=1 


Proof. First we show that £, o is well-dehned. Note that Hi vanishes at v = 0 by (SI) So, 
Remark [2~4] yields that d u H 2 (uo, lylv)) is exponentially localized. Combining the latter fact with 
Theorem |6.8[ -3 implies that the integral <f(A) converges and thus £ v o is well-dehned. It is well- 
known 02l Lemma 2.1.4] that, when the coefficient matrix depends analytically on a parameter, 
then the evolution is analytic in this parameter too. Combining this with Theorem 6.8 1,2 yields 
the first property. 
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A sketch of the proof of the second property 

Proving the second property is quite laborious, but not conceptually difficult. Therefore, we choose 
to give a sketch of the proof. The full proof can be found in the Appendix [C] It is sufficient to es¬ 
tablish that, for y e S , £ v o(d, y) is bounded away from 0 for A > 0 is sufficiently large. First, we 
prove that Q(A) remains bounded as A —* oo. Second, we show that a rescaled version of (3.7 1 has 
an exponential dichotomy on R for A > 0 sufficiently large. This enables us to construct suitable 
bases of the stable and unstable subspaces at 0, which will form the column vectors of an invert¬ 
ible matrix Tl{A ). Eventually, we estimate the product £ i; o(d, y) dct (77(d)) by a d-independent 
non-zero determinant for large A > 0. 


The third property 

Assume do is a simple zero of Gy o- Theorem 6.8 4 allows us to split off the singular part of Q(A) 
at d 0 . Indeed, we have for d e By 


Q(A) = 


d-d< 


'Pf + @ a (X), 


with tp, iJj and Q a (A) as in (6.6). Using the multi-linearity of the determinant, we expand for d e By 
and yeS 1 


<5.v,o(d, y) = det 


U a (A,y) + 


1 / 0 0 

A- A 0 \ pf 0 

2m 


T S (2L 0 ,0,A) 


= det(77, 


'a(A, 7)) + £ (fuflU do)) (p T [C„(do, 7) J m+1 ) ■ 


This concludes the proof of the third property. 


Remark 6.12 {Appearance ofy in the singular part). In the case m = 1, Proposition |4.6| shows 
that y appears as a factor in the singular part of the Laurent expansion of £«,o(-,y) at a zero d 0 of 
Gy o- Therefore, £ s ,o(-, y) has a pole at d 0 for some y e S 1 if and only if G s , oO, y) has a pole at A,, 
for all y e S *. Asa consequence, it is sufficient to check condition 2 in Proposition |3. 171 only for 
some y e S 1 . 


However, in the general setting of Proposition |6.1 1 -3 the principal part of the Laurent expan¬ 
sion of £i,o(-, y) is polynomial in y. So, it could happen that G St o(-, y) has a pole at d<> for all but a 
discrete set of y € S 1 . We expect that such a (non-generic) situation occurs precisely when d 0 is 
a limit point of the slow spectrum IJ rsS i ATfi.v.oF, y)). In this case, it is again sufficient to check 
condition 2 in Proposition [347] only for some y € S 1 . ■ 


It is possible to obtain expressions for the singular part of the Laurent series of £^,o(’, 7) 4 a zero 
of Gffl of higher multiplicity. Again, for simplicity of exposition, we treat only the case where 
d v G{uo, v, 0) is symmetric for all v € V. 


Proposition 6.13. Suppose d v G(uo,v,0) is symmetric for all v e V. Let d<> be a zero of Gyp of 
multiplicity ky. Letpyy = (vyy,qyy),i = 1,... ,ky,By and Xyj n as in Proposition 6.10 Define 
for i = 1 , ,ky and A € By 


Pi := 


-f 


d v H 2 (uo,v h (z))v U y(z)dz e C m , i/q := 


X oo 

w. 

oo 


y(z)*d u G{u 0 ,v h {z),0)dz € Mati xm (C), 
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and 


Qt 


X oo 

[d u H 2 (u 0 ,v h (z)) + d v H 2 (uo, Vh(z)) f VA„jn(z, T)] dz e Mat mXm (C), 

oo 


where r V,i j n denotes the upper-left (n x m)-block of the analytic (2n x 2 m)-matrix X,y xn . Let 
(ufx. A), pfx. A)), i = 1,... ,2m /;<? a fundamental set of solutions to the slow limit system 0 
For each subset re j I,..., k.i ( and each injective map cr : r —> { 1,..., 2777) denote by C TJT (A, y) 
the matrix obtained by replacing each cr{j)-tli column of the matrix 1L a (A, y) in (6.7) by (0, ipj)for 
all j € r. For all y € S the principle part of the Laurent series ofS s ^(-,y) at A 0 is given by 


Zttttv Z Z 


7=1 ^ rc{l,...,/t /U ), cr: 

#T=i 


J J 4*jM(r(j)(2'L 0» Ao) 


t—7=1 
injective 


det(C T)0 -(/l«,y)). 


Proof. This is analogous to the proof of Proposition [6JT} 3 using Proposition |6. 10|instead. □ 


6.4 An exponential dichotomy capturing the fast dynamics 

As mentioned in Section |6.1[ our goal is to apply the Riccati transformation on the linearized 
equations ( |3.3| ) in order to factorize £ £ into a fast and a slow part as in ( |1.2| ). However, according 
to Theorem |5.1| the Riccati transformation is only legitimate, when system ( |6. 1 [ ) has an exponential 
dichotomy on R. We can show that this is the case, whenever A is away from the zero set M(F>fp). 
Therefore, we introduce the following notation, before stating the result. 


Notation 6.14. Let A € (-Aq, 0) with Aq > 0 as in Lemma 6.3 For 5 > 0 we denote 


:= £a,o\ U B(A,8). 

,1eJV(£/,o) 


Theorem 6.15. Let A € (-Ao, 0) with Aq > 0 as in Lemma 6.3 Take 8 > 0. For s > 0 sufficiently 
small, system (6.1) lias for all A 6 an exponential dichotomy on R with constants Kf,py > 0, 
which are independent of s and A. 


Proof. Our approach is as follows. First, we show that system ( |6.1| ) has exponential dichotomies 
on the intervals [0,2L e ] and [-2L e ,0]. To establish an exponential dichotomy of (6.1) on an 
interval [a, 2 L s - a] for a > 0, we regard (|6.1|) as a perturbation of 


D 2 d x v = q 


d x q = (d v G(u p Jx). 0, e) + d) v 


v, q £ 


( 6 . 8 ) 


and apply Proposition |B.8| By choosing a > 0 independent of s, we can extend the exponential 
dichotomy to [0,2L £ ] by the extension Lemma|B3] Subsequently, we calculate the minimal open¬ 


ing between the stable and unstable subspaces at 0 of the two exponential dichotomies of (6.1) by 


comparing system (6.1) with the fast homogeneous limit system (3.5). We will show that, when¬ 
ever A is contained in X A ( >, this minimal opening is substantial. Therefore, an application of the 
pasting Lemma [BJ] gives an exponential dichotomy of (6.1) on [-2L £ , 2L £ ]. Finally, we apply the 
periodic extension Lemma B.6 to yield the result. 


An exponential dichotomy for system ( | 6 .8\ on R 

We show that the requirements in Proposition |B.4| are satisfied. First, by Lcmma[6R]thc coefficient 
matrix FI/, e (x, A) of system (|6.8[) is hyperbolic on RxCa, with eigenvalues bounded away from the 
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imaginary axis by some constant, which is independent of e, x and 4. Second, Remark [23] ini plies 
that is uniformly bounded onlx I ; \,o and that we have IflR/Ax, A) = 0( s[e) uniformly 
on R x Sa,o- So, by Proposition |B.4| system ( |6.8[ ) admits, provided e > 0 is sufficiently small, for 
A € Ea,o an exponential dichotomy on R with constants Kf r ,pf r > 0, independent of e and A. 

Exponential dichotomies for system f |6. /[ ) on [0, 2L S ] and on [ —2L e , 0] 

We want to apply Proposition |B.8| to system ( |6.8[ ). By assumption |(S3)| -3 there exists K x ,p\ > 0, 
independent of e, such that 

||<9 v G(« p , £ (x), v p>e (x), e) - d v G(u p ,s(x), 0, e)|| < K x e~^ mmU2L - 4 , (6.9) 

for x € [0, 2L e \. Now take xo > 0, independent of s, such that 

K x e~ MlX0 < 

36 K fr 

Estimate ( |6.9[ ) yields by roughness (Proposition |B.8[ ) that system ( |6. 1 [ ) has for every A e X A ,o an 
exponential dichotomy on [xo, 2 L s - xo] with constants independent of e and A. By the extension 
Lemma |B31 we can extend the exponential dichotomy of ( |6. 1 [ ) on [xo,2L e - xo] to [0,2L e ] with 
constants independent of e and A. An analogous treatment provides an exponential dichotomy of 
( |6.1| ) on \-2L g , 0] with constants independent of e and A. We conclude that ( |6.1[ ) has exponential 
dichotomies on both [0,2L e ] and [-2L e ,0] for every A £ I,\.o with constants Kfo,pfo > 0, inde¬ 
pendent of e and A. 


Comparison of system (6.1 ) with the homogeneous fast limit system (3.5) 


Take A e £a,o- Following Remark 2.5 we have for x e [s p , e p ] 


||d v G(«p je (x),v p , e (x),£) - 5 v G(uo,Vh(x),0)|| = 0(pf if ). 


( 6 . 10 ) 


Denote by T r (x,y,A) and Tp e (x,y,A ) the evolution operators of (|3.5|) and (6.1 1 , respectively. By 


Remark 2.3 the coefficient matrix of (|6.1[) can be bounded on R by some constant M > 0, indepen 


dent of e and A. So, system (6.1) has bounded growth with constants K x = I and p x = M. Now 


choosey = (3f/{AM) > 0. By Lemma B.2 and (6.10> we estimate, provided e > 0 is sufficiently 
small, for all x,y e [^log(e), -^log(e)] 


\\T r (x,y,A)-T f Jx,y,A)\\ < 1. 


( 6 . 11 ) 


The minimal opening for system (3.5) 

We recall some facts from Theorem 6.4 and Corollary |6.5| First, system ( |3.5| ) admits for A e Ca 
exponential dichotomies on both half-lines with constants K r ,p r > 0, independent of A. Second, 
the corresponding projections P r ±(x, A) are analytic in A. Third, the stable and unstable subspaces 
E ( ' + (0, /l) := P,+(0. /l)[C 2 "] and £“_(0, A) := kcr(£,_(0, A)) are complementary on L ; \,o. Therefore, 
Proposition B.14| implies that the continuous map p r '■ Ca —> [0, oo) given by the minimal opening 
rj r {A) = //(£,' + (0, A), £“_(0, A)) is bounded away from 0 on L ; \,5 by some constant ca ,<5 > 0. 


The minimal opening for system (6.1) 

Take A £ La.o- Denote by Pf± iE (x, A) the projections corresponding to the exponential dichotomies 
of (6.1 1 on [0,2L e ] and on [-2L £ ,0], respectively. By combining estimate ( 6.1 1 [ ) with Lemma 
there exists for each w e Ei, „(0,4) := Pf+A 0,4)rC 2 "l an element v £ E s r+ ( 0,4) such that 


B.ll 


P f+s (0,A):=P f+ A0,^[C Zn l 
11v - w\\ < (1 + K r ^)Kf 0 ^f°\\w\\. 


( 6 - 12 ) 
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Similarly, there exists for each w e E'j f (0, A) := kcr('Py_ E (0, /l)) a vector v e E“_(0, /l) such that 
(6.12) holds true. Therefore, Proposition |B. 14| -3 yields the estimate for A e Za ,<5 

o < c A ,a < Tj r (A) < q(E s f+E (0, A), E“_ £ (0, A)) + 4(1+ K r ^)K f0 ^°. (6.13) 


Application of the Pasting Lemma 


By estimate (6.131 and Proposition B. 14 -2 one deduces that, for e > 0 sufficiently small, E s , (0, A) 


and £y_ £ (0, /t) are complementary on Ha, a- So, the projection Pp E (A) onto E s f+e (0,A) along 


/+,£ V 


E u f _ (0,4) is well-defined on "Las- Moreover, by Proposition B.14-1 and (6.13) we obtain, pro- 
vided s > 0 is sufficiently small, the bound, 

ll^’/.eO'OII ^ n rat prc — > 

Tj(E s f+E (0, 4), E u f _ e ( 0, A)) c a, a 

for /I e Ha, a- Now, by the Pasting Lemma B.7 equation ( |6.1[ ) has an exponential dichotomy 
on [-2L £ , 2L e ] with constants depending only on CA,s,PfO and Kjq. Subsequently, by the periodic 
extension Lemnia [B~6| system ( |6.1[ ) admits an exponential dichotomy on R with constants Kf,pf > 
0 depending only on ca,s,P/o, E/o and M. Recall that M > 0 is a constant, independent of e and 
A, which bounds the coefficient matrix of ( |6. 1 [ ) on R x Ha,o- We conclude that pf and Kj are 
independent of s and A. □ 


For later purposes we treat the following corollary of Theorem 6.15 


Corollary 6.16. Let A € (—Ao, 0) with Ao > 0 as in Lemma 6.3 Take 5 > 0. For A e Z.ya, there 
exists a unique bounded solution x Pj n e (x, A) to the inhomogeneous matrix equation, 


dfV = JkiA*’ ^ + 3l 2 iAx), T e Mat 2nx2m (C), 
satisfying for all x 6 [ ~L r , L s ] 

\\^ in Jx,A)-X in (x,A)\\ =O(£?0, 


(6.14) 


uniformly in A € Ha, 6- Here, (3f > 0 is as in Remark 2.5 and Xi n (x,A) is the unique solution to 
established in Theorem 16. SI 


Proof. Let A e Ha, a- First, system (6.1 ) has by Theorem |6.15| an exponential dichotomy on R with 
constants Kf,py > 0, independent of e and A. Second, by Remark [23] we have for x € [ -s~ p , s~ p ] 

\\d v G(u v , s (x), v Pj£ (x), s) - d v G(u 0 , v h (x), 0)||, ||<9„G(n P;e (x), v P;£ (x), e) - d u G(u 0 , v h (x), 0)|| = 0(ePf). 


Now we apply Proposition |B. 12 to the inhomo geneo us equations (|3.6[) and (6.14): there exists a 


unique bounded solution '¥i n s (x,A) of equation (6.14) satisfying for all x e [-^e p ] 

\\^ in ,s(x,A)-X in (x,A)\\ =0(ePf), 


(6.15) 


uniformly in A e Ha, a- Flere, we have used Kj and pf do not depend on e and A and Xi n (-,A) is 
bounded on R uniformly in A e HA,a by Theorem 6.8 -3. Note that G vanishes at v = 0 by |(S1)| 
Hence, by |(S3)| -3 there exists K\,p\ >0, independent of e, such that for x € [-L s , L E ] it holds 

\\d u G(u p ^(x), v Pi£ (x), e)|| < (6.16) 


Combing estimate (6.16) with Proposition |B. 12| implies that there exists AT,q 2 > 0, independent 
of e and A, such that 

II^MII <K 2 e^ M , (6.17) 

for all x e [-L £ , L s \. Estimate ( 6.1 7[ ) together with the estimate in Theorem 6.8 -3 show that (6.15) 
actually holds for all x 6 [ -L e , L E \ uniformly in A e Ha, a- □ 
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6.5 Factorization of the Evans function via the Riccati transform 

As mentioned in Section [ 6 Tj an application of the Riccati transform enables us to reduce our linear 
stability problem ( |3.3| ) into two subproblems. This yields the factorization of the Evans function 
£ e into two factors £ Vjf; and £/. e , which can be linked to £/ ; o and £ v- o, respectively. The Riccati 
transformation is employed in the following theorem. 


Theorem 6.17. Let A e (-Ao, 0) with Ao > 0 as in Lemma 6.3 Take 6 > 0. For s > 0 sufficiently 
small, there exists a function U E : R x Ea ,<5 —> Mat 2 m x 2 «( ( C) such that for all A e Ea ,<5 and ye C 
we have the factorization, 


S e (A,y) = £,JA,y)£ f JA,y), 

with S s _ E ,&f E : Za ,<5 xC -> C given by 

£,, e (T, y) := dci(T sd J0, -L s , A) - yT sd J0, L E , A)), 
£ f ,s(A, y) := det(T/d, e (0, -L e , A) - yT, d J(). L s , A)). 

Here, Tsd.tAx, y, A) is the evolution of system, 

d xX = Vi(^n, £ (v, A) + Jl l2 , s (x)U E (x, A)) X , X € C 2m , 

and T f d jAx, y, A) is the evolution of system, 

d x a> - IJI 22 ,s(x. A) - sfsU E (x, A)Jln,e (x))to, to € C 2 ". 


(6.18) 


(6.19) 


Let /3f,p > 0 as in Remark 2.5 U E has the following properties: 

1. U E (-, A) is bounded uniformly in e and A on R; 

2. iff-. A) is 2L E -periodic for each A e Ea, 5 ; 

3. For all x € [-L e , L E ] it holds uniformly in A e 

\\Us(x,A) - Xi n (x,A)\\ = 0(sP f ). 


Here, M m (x, A) is the unique solution to (3.6 1 obtained in Theorem 6.8 


4. For all x e [-L E , L E ]\I E with I E := [-e p ,e p ] we have \\U F (x, A)|| = 0(s 3 ), uniformly in 
A e Xa j( 5 . 


Proof. Fix A € L t \ d . System (3.3 1 is clearly of the form (5.1 1 with coefficient matrices that are 


uniformly bounded in e on R by Remark |2.3| Furthermore, by Theorem |6.15 system ( 6 .1 1 has 
an exponential dichotomy on R with constants Kf,py > 0, independent of e and A. Hence, we 
can apply the Riccati transform from Theorem 


5.1 


and U E (x, A) £ Mat 2 mX 2 »(Q be as in Theorem 


5.1 


to (3.3 1 . Fet H e (x,A) 6 Mat 2 ( m +;i)x 2 ^ m + n )(C) 


<p(x) = H e (x, A)\Jj(x) transforms (3.3) into the diagonal system. 


By Theorem |5.1 \2 the change of variables 


d x f - 


yfs(Jln,s(x, A) + Jl 12 , £ (x)U E {x, A)) 


0 


Jli 2 , E (x, A) - xf£U E (x,A)Jln,e(x) 


f, ( 6 . 20 ) 


with f € C 2< '" + '' ) . The evolution T d , s {x,y, A) of system (6.20) is a block diagonal matrix with 
consecutively T s d, e {x,y, A) and Tf d , E (x, y, d) on the diagonal. Furthermore, H E {-,A) and U e {-,A) 
are 2L £ -periodic by Theorem|5.1|-5. Finally, as a product of two triangular matrices with only ones 
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on the diagonal, the determinant of H e {x, A) equals 1 for every x e K. This enables us to factorize 
the Evans functions S E as follows 

£ e (A, y) = det (h e { 0, A) [T d , e ( 0, -L E , A) - yT d J0, L s , A)] H E {L E , d)" 1 ) = £ s , e (A, y)£ f JA, y). 


Properties of U e 

The first two properties are immediate by Theorem |5.1[ -1 and Theorem |5.1| -5. Let ffV Jj£ (x, A) be 
the unique solution to (6.14 >. By Theorem |5.1| -3 and Corollary |6.16| it holds for a e [-L £ ,L £ ] 

II U E (x, A) - X in {x, d)|| < II U E {X, A) - TV £ (x, .1)11 + X) ~ Xin(x, Till 

= 0(e min|1/4 ^/ 1 ), 

uniformly in d € La, a- Here, we have used that Kf and py arc independent of e and d. This 
settles the third property, since we may without loss of generality assume /?/ < 1/4. For the fourth 


property we use the method of successive approximation. Note that G vanishes at v = 0 by (SI) 
Hence, by |(S3)| -3 there exists K\,p\ > 0, independent of e, such that ||J?l 2 i,£(x)ll < K\e~^'^ for 
x € [-L £ ,L £ ], Using the latter two lines, we approximate U e {x,A) successively for three times 
with Theorem |5.l| -4. This yields \\U e {x,A)\\ = 0(e 3 ) for x e [~L e , L E ]\I E uniformly in d e Ea,< 5 - 
Here, we have used again that Kf and pj are independent of e and A. □ 

6.6 The factors of the Evans function and the three singular limit problems 


We derive from Theorem 6.17 that our linear stability problem (3.31 diagonalizes into two sub¬ 


problems (6.181 and (6.19). This diagonalization yields the splitting of the Evans function £ e into 
two factors £ se and S/ E . By relating (6.18) to the slow limit problem (3.7) and (6.19) to the ho¬ 
mogeneous fast limit problem (3.5 l, we can link £ S E to the slow reduced Evans function £ A p and 
&j E to the fast reduced Evans function fi/Q. 


We start with the relation between (|3.7[) and (|6.18 1. Note that the transformation matrix U, 


is 


related to the fast inhomogeneous limit problem (3.6) by Theorem 6.17 -3. This is how problem 


merges into the slow reduced Evans function o via system (6.18). This is reflected in the 


proof of the following result, which links S s ^ E to the slow reduced Evans function £ v .o- 


Lemma 6.18. Let A e (-Aq,0) with Aq > 0 as in Lemma 6.3 Take 6 > 0. Let <5 (£ be as in 


Theorem 6.17 For e > 0 sufficiently small, there exists p s > 0, independent of e and A, such that 


we have the following estimate, 


£,UXy) = £sfi(A,y) + 0{^), 


uniformly in A 6 Ea ,<5 and y e S l . 

Proof. Let d e Ea ,<5 and y e .S' 1 . Our approach is as follows. We split the coefficient matrix in 


system (6.18) corresponding to their decay behavior outside the pulse region, i.e. we write 

JlnAx, A) + WnAxWsix, A) = Si, £ (x, d) + S 2 . E (x, A), (6.21) 

with 

0 


S 1;£ (x,d):= 
S 2 , e (x,d) := 


D T 1 

s (d u H\ {tip A x ), 0, £) + A) 0 

0 0 


+ ^12 Ax)U e {x,A), 


BiAx) 0 

BiAx) := (fHiJupJx), VpAx)) + £ (d u H](iipJx), v Pi£ (a), e) - d n H \(itpjx), 0, e)) . 
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Note that A) and B 2e (-, A) are bounded on IK. uniformly in e and A e £ A S by Theorem 6.17-1 


and Remark 2.3 The splitting gives rise to an intermediate system, 

d x x = y[sB\, e {x,A)x, x e 


■^2 m 


( 6 . 22 ) 


The intermediate system (6.22 > helps us to draw the connection between system (6.1 8|) and the 


slow limit system (3.7). On the one hand, we will show that (6.1 8|) and (6.22) are closely related 


by using variation of constants. On the other hand, by performing a suitable coordinate change the 


evolutions of systems (3.7) and (6.22 1 are close to each other. This enables us to approximate <5 V 
with £ s? o on £ Aj5 . 


Some estimates 

Denote by Ti SjE (x,y, A) the evolution of system (6.22). Let p > 0 be as in Remark [23] By Lemma 
B.2 we have for each x,y e I E - \ -s~ p , s~ p \ 


I \T s d,s{x,y, A) - 7||, || T is , E (x,y, A) - 7|| = 0(s 1/2 ~P), 


(6.23) 


uniformly in A e Next, we will pay attention to the term B 2 , E in the expansion (6.21 1 . 

First, note that H 2 vanishes at v = 0 by |(S1)| Hence, assumption |(S3)| -3 implies that there exists 
K\,p.\ > 0, independent of s, such that ||7?2, e (x)|| < for x e [-L E ,L e \. Putting the latter 

and Theorem 6.17|-4 together, we approximate for v € [-L e , L E ]\I E 


\\B2 Ax,A)\\ = 0(s 3 ), 


(6.24) 


uniformly in A e On the other hand, one readily establishes that (6.22) has bounded growth 
on [- L e ,L e ] with constants K] £ ~ l/2 , sp\ > 0, where K\,p\ > 0 are independent of A and s. 


Therefore, another application of Lemma B.2 gives for x,y e [0, L r \ \ /,. 


\\T sd , s (x,y,A)-T is , E (x,y,A)\\ < K,s~ xl2 L E 


sup Vi||S 2 , £ M|| e p ^=0(s 2 ), 

xE.[0,L £ ^\jl e 


(6.25) 


uniformly in A e Z A ^. 


Connection between (6.18 1 and (6.22): a variation of constants approach 
The variation of constants formula gives 

Ts d J(L L e , A) = T is , E ( 0, L e , A) - s[e f T is J0, z, A)B 2 Jz, A)T sd , E (z, L e , A)dz. 

Jo 


Using the bounded growth of (6.22), we have by estimates (6.24) and (6.25) 

r eP 

T sd , e ( 0, L e , A) = T is J0. L e , A) - \fr I T IS J0. z, A)B 2 Jz, A)T fd Jz, L e , A)dz. + 0(s ) 

Jo 

= r + , E (J)T is , E (0,L E ,A) + O(s), 

uniformly in A e £ A ,^, with 


(6.26) 


T + , e (A) ■= 1 


Vi f 

Jo 


TtsA 0, z, A)B 2 . e (z, A)T sd , £ (z, s~ p , A)dzT l sAs~ , \ 0, A). 


Using (6.231, we derive 


ps p 

T + ,M) = / - Vi I 'B 2 Jz, A)dz + 0(s l ~ 2p ) 
Jo 


(6.27) 
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uniformly in A e Ea,< 5 - Write the (2nx2m)-matrix X m (x, A) as a composition of four block matrices, 
where < V,„(x, 2) is the upper-left n x m-block. Using Remark [23] Theorem |6.8[- 3 and Theorem 


6.17 -3, we approximate (|6.27|) as 

T + JA) = I + ^ f 
Jo 


0 0 

d u H 2 (u 0 ,v h (x)) 0 

0 


+ 


0\( f V in (x,A) 0 
* 0 


dx + 0(e l,2+fI f (6.28) 

+ 0(E l/2+fJs ), 


d v H 2 {u 0 ,v h (x )) 0 

I 0 

Ve Jo" [d„H 2 (uo, Vh(x)) + d v H 2 (u 0 , v h (x))*V ( > l (x, d)] dx I 

uniformly in A e Z,\, ( >, where p s = (if - p. Without loss of generality we may assume 0 < p < /3f 
so that p s > 0. 


Connection between system (3.7) and (6.22) 

We apply two operations on system (6.22 >. First, we perform the coordinate change = C E x, 
where C E := (q € Mat 2 mX 2 m(Q- Second, we switch to the large spatial scale x = ex. The 
evolution of the resulting system from 0 to L e can be approximated by 7^(0, Lo, A) by combining 
Lcmma[B.2| with (S3) 1,3. Thus, we obtain the following estimate, 


II C^T^CO, L e , d)C e - T s ( 0, L 0 , d)|| = <9( >£), 


(6.29) 


uniformly in A e 


Comparing S S E to the slow reduced Evans function S s p 
Plugging (6.28 1 and (6.291 into ( |6.26 ) yields 

TsdA 0, L e , A) = C E [t+(T)7T(0, L e , A) + 0(^ )\ C~\ 

uniformly in A e with 


(6.30) 


T+(d) := 

Similarly, we derive 


/ 0 

- f 0 °° [d u H 2 (uo, v h (x)) + d v H 2 {u Q , Vh(x)) f Vj n (x, d)] dx I 


T sd ,s( 0, -L e , A) = C E [t_(T)7T(0, -L e , A) + C~\ (6.31) 

uniformly in A e with 

( I 0 

T ^ ' v f°oc [d»tf 2 (wo, VhW) + d v H 2 (u (h v h (x)) < V i „(x, d)] dx I 

First, by Liouville’s Theorem we have det(7~( s (0, To, A)) = I. Second, the identity Y+(d)~ ! Y_(d) = 
Y(/l) holds true, where T(A) is defined in (3.8 1 . Putting these two items together, we estimate Ca¬ 
using ( |6.30| ) and ( |6.3 1 [ ) 

£s, s (A, y) - det(7~^ e (0, -L e , A) - YT sd , s (0, L e , A)) 

= det (r-(A)T s (2L 0 , U, A) - 7 T + (A)T s (0, L 0 , A)) + O(rf ') 

= £,,o (A,y) + 0(tf‘), 

uniformly in A e £a ,s and y e S 1 . This estimate concludes the proof. □ 


42 





















It remains to link Sf E to the fast reduced Evans function £ 


'/■ o- 


Lemma 6.19. Let A e (-Ao,0) with Ao > 0 as in Lemma 6.3 Take 5 > 0. There exists p p > 0 
such that, for s > 0 sufficiently small, there is a map h E : Xa,) —> C satisfying 

0 < \h E (A)\ = 0(e-^), 

£f, s (A,y)h E (A) = (-y)"£/,o(d) + 0{sTr), 

uniformly in A € Xa ,<5 and y € S l . 


Proof. Our approach is as follows. First, we show that system (6.19» has an exponential di¬ 
chotomy on R, if A is in XaR ecall that, the fast reduced Evans function £/,o is defined in terms 
r (A) of the stable and unstable subspaces at 0 of the homogeneous fast limit system 


of bases B u s 


(3.5 i. By comparing system (6.191 to (3.51, we are able to construct bases Bf (A) of the (un)stable 
subspaces at 0 of (6.19», which are close to B l f s (A). By tracking the bases Bf '\A ) either forward 


or backward, we construct bases of the (un)stable subspaces at ±L e of (6.191. These bases at ±L E 
will form the column vectors of a matrix LL E (A), which connects £y e to £yo- 


An exponential dichotomy on R of system ( |6. 19 ) 

Let A e Xa ,<5 and y e .S' 1 . First, system (6.1 1 has by Theorem |6.15 an exponential dichotomy 
on R with constants Kf,py > 0, independent of s and A. Second, U s {x,A) is bounded on R uni¬ 
formly in e and A e Xa,,s by Theorem 6.17 1. Therefore, Proposition |B.9| yields that (6.19 1 has an 
exponential dichotomy on R with constants Kfd,p/d > 0, independent of e and A. Denote the cor¬ 
responding projections by Pfd,e(x, A), x e R. Since the coefficient matrix of (6.19) is 2L £ -periodic 


by Theorem 6.17 -2, the projections Pfd.e{-, A) are also 2L £ -periodic by |[9[ Proposition 8.4]. 


Comparing system {6.19) to the homogeneous fast limit system (3.5) 


Let M > 0 be a bound of the coefficient matrix ^ 122 , 0(^5 A) of equation ( |3.5[ ) on 
Remark |23] and Theorem 6.17 1 we have for a 6 f = [ -s~ p 


11^22 ,s(x,A)- sfeU e {x, A)5Rn,E{x) 


,£~ p ] 

Jh2,o(x,A)\\=0(s P0 ). 


X Xa,<5. By 


(6.32) 


with po := minj j3f, 1/2). Denote by T r (x,y,A) the evolution of (3.5) and take^ := po/(4M) > 0. 


By Lemma B.2 and estimate (6.32 1 it follows, provided that s > 0 is sufficiently small, for all 
x,y e \y logfg), -y logA;)] 

\\T r (x,y, A) - T fd , s (x, y,A )|| < 1. (6.33) 

By Corollary |6.5| system ( |3.5[ ) has for A e Xa.o exponential dichotomies on both half-lines with 
constants K r ,p r > 0, independent of A. Let Bf\A) be as in Theorem 6.4 Since Xa,o is bounded 
and B l f \-) is continuous, there exists a constant Ca > 0 such that ||fi/’ ,v (T)|| < Ca for A € X/ 


-A,0- 


Now, combine estimate (6.33) and Lemma [B . 11 1 there exists, for s > 0 sufficiently small, bases 
B l f s : Xa,) -» Mat 2nx „(C) o\P fd J(),A)\C 2n \ = ^[C"| and ker(^ /d , £ (0, A)) = B“[C n ], such that 


l|B“’ s (T) - R^(T)|| = Oca**), 


(6.34) 


uniformly in A e Xa,). Via (6.34 1 we establish, provided that s > 0 is sufficiently small, the bound 
ll# , (A)|| < Ca for A € X A , 5 . 


Comparing £/ £ with the fast reduced Evans function £yo 
Define for A e Xa,, 


TIM) ■.= (T Me (-L s , 0, A)B U JA), T fdtE (L s , 0, A)Bf(A)) . 
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(6.35) 


Sine ePfd, £ {', T) is 2L £ -periodic, the first n column vectors o\'H,(A) form a basis of ker(!P/^ £ (L £ , A)) 
and the last n column vectors form a basis of P f d ^ E (L £ , A)[C 2n ], Thus, H S (A) is invertible. By 
Hadamard’s inequality we have dcl('HJA)) = 0(e~ 2mfdLe ), uniformly in A £ L ; \,o- Moreover, 
since Pfd,s(-, A) is 2L e -periodic, we estimate 

HT/^fO, L s , X)T fd A-U, 0, T)B £ (T)|| < K 2 fd C A e~ 2 ^, 

\\T fd , E (0, -L E ,A)T fd , s (L £ ,0,A)B s E (A)\\ < K 2 fd C A e~ 2 ^, 

for A e S ; \.<5. We combine estimates ( |6.34 ) and ( |6.35 ) and derive 

(r fd , e (0, -L s , A) - yT fd J(l L £ , A)) HA A) = { B u r (A ), yB s r (A)) + O(N^), 

uniformly in A e Ea,<; and y e S 1 . Taking determinants in the latter matrix equation yields 

£/>(T, y) det(H s (A)) = (- y)"£/,o(T) + O(^), 

uniformly in A e Ea ,<5 and y £ .S' 1 . Thus, defining liJA) := dct(7£(T)) concludes the proof. □ 


Remark 6.20 (The connection between and Sffi). In the proof of Lemma 6.19 the connec¬ 
tion between &f s and £yo is given by the matrix H s . This idea is taken from the proof of l46l 
Theorem 2], However, the context in ||46|| is different. Here, one shows that the eigenvalues of a 
periodic boundary value problem are exponentially close to the eigenvalues of the corresponding 
unbounded problem. ■ 


6.7 Application of Rouche’s Theorem 


In contrast to the estimate achieved in Lemma 6.18 we need to rescale &f £ in Lemma 6.19 by 
an exponentially small quantity h E in order to relate it to the e-independent fast reduced Evans 
function S/o. This quantity prevents us from directly estimating the Evans function £ e by the 


reduced Evans function £q (A,y) — (-y)''£ v .o(d, y)£/.o(,1), using the estimates in Lemmas 6.18 


and 6.19 Nevertheless, it is still possible to compare the zero sets of S e and So using the classical 
symmetric version of Rouche’s Theorem due to Estermann. This will conclude the proof of our 
main results. 


Proof of Theorem 13.101 Let A € (—Ao, 0) with Ao > 0 as in Lemma [63] Take 6 > 0 sufficiently 
small such that T c By Proposition 16.1 111 the slow reduced Evans function <S v ,o(vy) is 


meromorphic in the interior of L and analytic and non-zero on T. Moreover, by Theorem 6.4 the 
fast reduced Evans function £yo is analytic on T and its interior and is non-zero on T. Hence, for 
e > 0 sufficiently small, we have for all A e T and y € S 1 


\S E (A,y)-S 0 (A,y)\ 

< (1 - h e (A))\S e (A,y)\ + | S s , £ (A,y)S f , £ (A,y)h £ (A) - (-y)"£ i ,o(d,y)£ / ,o(d)| 
<|£ e (d,y)| + |£ 0 (d,y)|, 


(6.36) 


where we have used Theorem 6.17 and Lemmas 6.18 and 6.19 The result follows by an applica¬ 
tion of the symmetric version of Rouche’s Theorem. □ 


Proof of Theorem l3.111 Let A e (-Ao, 0) with Ao > 0 as in Lemma 6.3 Let y e S 1 . By Propo- 
sition |6.lT -1,2 and Theorem 6.4 the map fi y : I,\.o —> C given by £ r (d) = (-y)''fi s .o(4, y)£/,o(d) is 
non-trivial and analytic. The zeros of £ y are given by the discrete set N y := N(Sj\o) U N(S s iA 'ey))- 
Now, take 5 y > 0 sufficiently small such that for every 4 0 € N y the disc B(d 0 ,d r ) contains no 
A e N y with A ± To. Take 0 < 8 < 8 y and T<> € N y . For all A £ d B(T<>,d) we have by construc¬ 


tion that £ r is non-zero. Hence, for e > 0 sufficiently small, inequality (6.361 holds true for all 
A £ 5R(T 0 , 8). The result follows by an application of Rouche’s Theorem. □ 
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Proof of Proposition [3TT6] Let A e (-Ao,0) with Ao > 0 as in Lemma 6.3 Let A 0 be a simple 
zero of Sffi. By hypothesis, there exists t?i > 0 such that 


fl(/L,5i) c {/l 0 } U 


Sa.o\ 


U A(£ s ,o(-,y))uA(£/,o) 


K~yeS 1 


Take 0 < 6 < 5\. The disc B(A 0 ,8) contains no zeros of <S v o(-, y) for every y e .S' 1 and precisely 


one simple zero A 0 of Hence, for e > 0 sufficiently small, inequality (6.361 holds true for all 
A 6 dB(A <> ,8) and y e .S' 1 . The result follows by combining Rouche’s Theorem and Proposition 

I6TT1 3. □ 


Remark 6.21 (Rate of convergence). The technical results Lemmas 6.18 and 6.19 seem to provide 
a rate at which the spectrum cr(_£ e ) converges to the spectrum in the singular limit. However, the 
approximations in these lemmas are only valid away from the zeros of the fast reduced Evans 
function <5yo '■ So, one can only deduce that spectrum converging to 


U Af(6 J>0 (-,y)) 

.yeS 1 


\ N(6 f , o), 


does this at an algebraic rate of order 0(£^ s ). By making the parameter 8 appearing in the proof 
of Lemma [6T9] dependent on e, it may be possible to derive an overall rate at which the spectrum 
cr(£. E ) converges to its singular limit spectrum. However, this is beyond the scope of this paper. ■ 


7 Concluding remarks 

7.1 Discussion 

As mentioned in the introduction, our factorization method via the Riccati transformation of the 
Evans function offers one unified analytic alternative to both the elephant trunk procedure de¬ 
veloped by Alexander, Gardner and Jones |[ll [22j] and the NLEP approach of IT3l fl2ll - that both 
have a geometric nature. It is worthwhile to compare and discuss the links between these methods. 

Moreover, the present work can be seen as the natural generalization of the spectral analysis 
f53l of periodic pulse solutions in the Gierer-Meinhardt equations to periodic pulse patterns in 
the general class of ‘slowly nonlinear’, (m + n)-component, singularly perturbed reaction-diffusion 
equations ( |2.3| ). Recently, a similar generalization to slowly nonlinear 2-component systems for 
homoclinic pulses has been developed in llT8l . In that sense, the present paper stands in the tra¬ 
dition of ||T2l[T8]|53l. However, our spectral analysis differs fundamentally from the analyses in 
these works, which rely eventually on the geometric approaches developed in HI |T3j [ 22l . 

7.1.1 Relation to the elephant trunk procedure 

Consider a localized pulse solution to a 2-component, singularly perturbed reaction-diffusion 
equation. When the linear stability problem (LSP) has a slow-fast structure, it is a general phe¬ 
nomenon that it decouples outside the pulse region due to exponential decay of the solution to 
the asymptotic background state. This yields a decomposition of the solution space into three 
subspaces V s± © V c± © V u± at both sides (+) of the pulse region. Here, V s± consists of fast ex¬ 
ponentially decaying solutions. Similarly, V u ± consists of fast exponentially increasing solutions. 
Lastly, V c± consists of solutions that evolve slowly. In the sense of ll39l . one could say (LSP) 
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admits exponential separations with respect to the decompositions © V c± © V u± . The difficulty 
is to ‘glue’ the subspaces V.+ and V._ for • = u,s,c together, yielding an exponential separation 
of (LSP) on the whole line. Eventually, this induces a factorization of the Evans function in a fast 
and slow component. 

Gardner and Jones achieved this in If22l by considering (LSP) in projective space. When (LSP) 
is asymptotically of constant coefficients type, one can first obtain stable and unstable bundles. 
Subsequently, these bundles are split into fast and slow (un)stable subbundles. The elephant trunk 
lemma is used to track the fast (un)stable bundle through the pulse region. By the control on the 
fast subbundle, it is possible to approximate the dynamics of the slow (un)stable subbundles. Even¬ 
tually, this yields a (1,2,1)-exponential separation of (LSP) on R. Note that the 2-dimensional 
center direction corresponds to the slow stable and unstable subbundles. In our stability analysis, 
the Riccati transformation plays the role of the elephant trunk lemma - see Section [63] This trans¬ 
formation yields an ( n , 2m, /tj-exponcntial separation on R of (LSP) as long as we arc not close to 
the eigenvalues of the fast singular limit problem. 

Although the proof of the elephant trunk lemma has been worked out in full detail for some 
specific 2-component models Ifl3l [T9l [22l 144 1 only, it is widely accepted that the method can be 
followed for a larger class of systems. However, there are some limitations. Lor instance, the 
elephant trunk lemma is only suitable for linear stability problems that have an asymptotically 
constant coefficient matrix. This is neither a restriction for slowly linear systems as the classical 
Gray-Scott and Gierer-Meinhardt models nor for homoclinic pulses on R. However, the linear 
stability problem associated with spatially periodic patterns in slowly nonlinear systems exhibits 
non-autonomous behavior in the background state on its domain of periodicity - and thus does 
not approach a constant coefficient matrix. This prohibits the application of the elephant trunk 
procedure. Moreover, the elephant trunk lemma is only capable of tracking the ‘most unstable’ 
fast solution, which corresponds to the (simple) eigenvalue of largest real part of the asymptotic 
coefficient matrix. Therefore, it is unclear - how to obtain the exponential separation with the ele¬ 
phant trunk method in the multi-dimensional setting n > 1. 

Lurthermore, there is a major difference in the mathematical framework used in [[0 [22] and our 
work. The framework in fH [22j has a highly geometrical character, whereas our method is of a 
more analytical nature. Alexander, Gardner and Jones track solutions via vector bundles formed 
from the projectivized (LSP). This has the advantage that the generated bundles have a clean and 
natural characterization as s tends to zero, whereas the actual solutions of (LSP) become singular. 
On the other hand, one could argue that exponential dichotomies provide a natural framework to 
capture the dynamics of (LSP) being a non-autonomous linear system, which depends analytically 
on the parameter A. The Riccati transformation is naturally formulated in terms of exponential 
dichotomies and is explicit in terms of the coefficient matrix of (LSP). Therefore, the exponen¬ 
tial separation of the solution space is much more explicit than in IHIZD, which shortens proofs. 
Linally, it is interesting to remark that in both the approach initiated by Alexander, Gardner and 
Jones and our method we need an a-priori e-independent estimate on the sector containing the 
spectrum. Our proof of this fact in Section [672] forms an analytical counterpart to the geometrical 
proof provided in fU Proposition 2.2] and If22l Lemma 3.3]. 

7.1.2 Relation to the NLEP approach 

Based on the geometric methods of Alexander, Gardner and Jones Hl[22]], the NLEP approach was 
developed in the context of the stability of homoclinic /V-pulsc patterns in the Gray-Scott equa¬ 
tion Ifl3l and Gierer-Meinhardt-type models Ifl2l . This method established the approximation of 
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the Evans function by the product < |1,3[ > of an analytic fast reduced Evans function and a mero- 
morphic slow reduced Evans function and provided explicit analytic expressions for both factors. 
The NLEP approach was extended to the spectral analysis of spatially periodic pulse patterns (in 
semi-strong interaction) in the generalized Gierer-Meinhardt equations in ll53l and to the stabil¬ 
ity of heteroclinic and homoclinic multi-front patterns in 2- and 3-component bistable systems 
of FitzHugh-Nagumo-type Ifl4l l55l . Moreover, the method has recently been generalized to the 
stability of homoclinic pulses in slowly nonlinear systems in Ifl8l 1551 . In each of these works, 
the fast and slow reduced Evans functions are interpreted geometrically in terms of fast and slow 
transmission functions that encode the passage of specially selected fast and slow basis functions 
over the fast pulse regions. The expressions for the slow transmission functions include Melnikov- 
type components. The meromorphic character of the slow reduced Evans function generates the 
zero-pole cancelation mechanism - also called NLEP paradox - in each of these models. The 
spectral analysis for periodic pulse solutions developed here shows that these phenomena occur in 
a broad class of multi-component singularly perturbed reaction-diffusion systems. 


Although the present work stands in the tradition of Ifl3l fl2l 14, 18:, [55, (53j, the methods dif¬ 
fer fundamentally. Unlike these works, our analysis is based on an intrinsically analytic reduction 
method. This has the advantage that our spectral analysis allows for non-autonomous behavior of 
the linear stability problem outside the pulse region - a crucial extension in the case of spatially 
periodic patterns in slowly nonlinear systems. This extended applicability of the present method 
also plays a role in the spectral analysis of homoclinic patterns: for instance a wv-term in the v- 
component of the n - m - 1 version of <o> cannot be allowed in Ifl8l . whereas our stability 
analysis can indeed handle such terms. Moreover, unlike in the present work, the singular limit 
problems appearing in lfT3lfl2l 11411181, 55 ]l53ll are scalar, which significantly simplifies the analy¬ 
sis of these problems. In |f j~3l f~2 531 the slow and fast reduced Evans functions can be explicitly 
computed in terms of hypergeometric functions, while in lfT4l l55ll the stability of the (multi-)fronts 
is determined by spectrum near the origin, so that the relevant reductions can be determined in a 
relatively straightforward manner. An extensive analysis of the multi-dimensional singular limit 
problems, as we did in Section 6.3 is thus not necessary in these cases. Finally, there also is an 
important similarity between the geometric methods in the literature and the analytic approach 
developed here. Although the methods are most often applied to patterns that exhibit a reversibil¬ 
ity symmetry, this symmetry is not essential for the application of the methods - see f55il for 
an example. Nevertheless, such a symmetry in general simplifies the analysis - as mentioned in 
Proposition |4.4| 


7.2 Future directions 


Perhaps the most pressing question is how the fine structure of the ‘small’ spectrum around zero 


can be determined (Section 3.4). This fine structure is crucial to decide upon spectral stabil¬ 
ity (Corollary 3. 17| ) and, eventually, nonlinear diffusive stability (Remark |3.19[ ). In Il47l the fine 
structure is presented for nearly homoclinic wave trains and in lf53l for periodic pulse solutions of 
semi-strong interaction type in the setting of the slowly linear Gierer-Meinhardt equation. How¬ 
ever, the methods in these two papers are not directly applicable to our situation - see Remark 


3.18 Nevertheless, preliminary investigations indicate that the analytic methods of PiTll may be 
extended to the present semi-strong interaction regime by separating fast from slow dynamics via 
Riccati transformations. This is the subject of work in progress. 


Another direction is the nature of destabilization of spatially periodic pulse solutions. Recent 
research ifTTl |54j shows that destabilization mechanisms can be rather complex when periodic 
patterns approach a homoclinic limit. While increasing the wavelength, the character of the desta- 
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bilization alternates between two kinds of Hopf bifurcations. This phenomenon is called the ‘Hopf 
dance’ lfT71[54l . It has been analytically established in (slowly linear) Gierer-Meinhardt models 
in m and recovered by numerical methods in the generalized Klausmeier-Gray-Scott model 
mm. The latter observations suggest that it is a persistent mechanism - at least in slowly linear 
models - especially since the systems considered in Il54l include nonlinear diffusion terms. How¬ 
ever, the analysis in ifTTll suggests that slowly nonlinear dynamics of the models may a priori have 
a decisive impact on the appearance of the Hopf dance. Both the Hopf dance as well as the ‘belly 
dance’ IfTTll - an associated higher order phenomenon - can be analyzed in the general slowly 
nonlinear setting of ( |1.1| ) by the methods developed here - see iflOl . 

Finally, we note that the analysis developed in this work may be used as the foundation for an 
analytic study of the bifurcations exhibited by spatially periodic patterns in semi-strong interac¬ 
tion - patterns that thus are ‘far from equilibrium’. In other words, the developed explicit insights 
in the linear stability is a key to understanding the weakly nonlinear dynamics of patterns as they 
have become spectrally unstable. A first - and fundamental - step in this direction has been taken 
in 115711 . in which a normal form approach associated with a Hopf destabilization of homoclinic 
pulses in the n - m - 1 version of <o> is developed. Unlike known classical slowly linear ex¬ 
amples such as the Gray-Scott and Gierer-Meinhardt models, the Hopf bifurcation for homoclinic 
pulses can be supercritical. It can even be the first step in a sequence of further bifurcations that 
leads to complex (amplitude) dynamics of a standing solitary pulse - as observed in the simulations 
of l58l . It should be noted that the subcriticality of the Hopf bifurcation in Gierer-Meinhardt-type 
models was a ‘conjecture’ based on numerical evidence until the work in 1571 . Both the above 
described Hopf dance and the fact that the pulses that together form the spatially periodic pat¬ 
terns are in semi-strong interaction, indicate that the weakly nonlinear dynamics of these patterns 
beyond their destabilization may be very rich. 


A Proof of existence of stationary, spatially periodic pulse solutions 


Proof of Theorem 12.1 II We adopt the notation of assumptions |(E1)| and |(E2)| Recall that (12.41) 


is /(-reversible by Remark 2.2 Our proof is based on the fact that every orbit that crosses the 


space ker(7 - R) twice, must be a closed loop. Therefore, we start with a ‘good’ set of initial 
conditions Z c ker(7 - R) and track these conditions under the forward flow of \2A\ with the 
aid of an appropriate Exchange Lemma. We will show that the tracked trajectories remain close 
to the singular orbit consisting of the segments i// s (x), x e [0, 2 Lq] and <p\fx, no), x e R, where f h 


is defined in Remark 2.9 In particular, we establish that the union of trajectories starting in Z 
intersects ker(7 - R ) transversally in some point P pE , which lies close to 0h(O, no). Finally, the 
desired periodic solution is the one that starts in P pE . 


A good set of initial conditions 

Denote by e,-, i e 1,..., m the unit basis of R m . Let Id be 
vectors e \,..., e/„_i, e^+i,... ,e m , where i* is as in 
manifold, 


(E2) 


the (m x (in - l))-matrix with column 
Consider the (m + n - l)-dimensional 


Z ■= {(Ml + n, 0, V, 0) : u e T/[R m_1 ], vefjc ker(7 - R). 


The intersection of Z and A1 corresponds to, 

P 0 := {(ni + u, 0) : u e '77[R m “ 1 ]} c ker(7 - R s ). 
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By assumption (E2) Po becomes under the forward flow of the slow reduced system (2.6 1 an 
m -dimensional manifold PI, which intersects TL transversely at (mq, -ff(iio))- Indeed, we have 


0 + det 0,(0, U) 


= det 


~d u J (mo) 

I -J{u o) 

d u J(u 0 ) 77i(w 0 ,0) 


0 PI 
H\(u\, 0) 0 

PI 


®s(L 0 ,0) 


0 


We have used here that O s (Lo, 0) induces an isomorphism between the tangent spaces of P* Q at 
(i/i,0) and at (wo, —JT(u o)) and that the determinant of O s (Lo,0) equals 1 by Liouville’s Theorem. 


Putting system (2.4) in Fenichel normal form 

Let Alo be a compact 2/w-dimensional submanifold of At so large that Alo serves as a neighbor¬ 
hood of il/ s (x), x € [0, Lo] an d of the projection (wo, C' Vh(z, uf))dz), x e R of (p\fx, uq) on 

At 


By assumption |(S2) Alo is normally hyperbolic. So, according to Fenichel’s theory li2Tl Theo¬ 
rem 9.1], Alo perturbs, for s > 0 sufficiently small, to a manifold Al £ , which is diffeomorphic 


to Alo and locally invariant for the dynamics of (2.4 1 . Since Alo is itself locally invariant for the 
dynamics of (2.4), one readily establishes that there exists a constant Co > 0 such that Al e has 


Hausdorff distance 0(e c °/ £ ) from Alo - see also IfTTl Theorem 2.1]. 

By OT1 Proposition 1] there exists a C 1 coordinate change ( u,p,v,q) —» ( a,b,c ), which brings 


system (2.4 1 into ‘Fenichel normal form’, 


d x a = A(a, b, c, e)a 
d x b = B(a, b, c, e)b 
d x c = sK{c,e) + C(a,b, c, e)(a® b) 


a,b e R n ,c e 


»2 m 


(A. 1) 


in an e-independent neighborhood T) of Alo, where A, B, K and C are continuous, A is a vector 
A,B are square matrices of order n and C is a tensor of appropriate rank. In particular, 


in 


i>2 m 


C(a 0 b) is bilinear in a and b. In ‘Fenichel coordinates’ Al e is given by a = b = 0 and the local 
stable and unstable manifolds W l f\M E ) HD of Al e are the spaces b = 0 and a = 0, respectively. 
Since R maps W"(Al £ ) to VF,!( Al f: ), ker(7 - R) n T) is contained in the space a = b. Finally, system 


d x c = K{c, 0) corresponds to the slow reduced system (2.6 1 . 


Application of the Exchange Lemma 

By the latter paragraph Z. c ker(7 - R) intersects the local stable manifold WC(Alo) n D of the 


fast reduced system (2.5) transversally at (wi, 0,0,0). Moreover, the reduced slow flow (2.6) on 


At is not tangent to Pq at (wi,0). We conclude that the conditions for the Exchange Lemma 
Section 2.5] arc satisfied. 


Denote by P s c Al £ the (m - l)-dimensional manifold, where Z and the local stable manifold 
W s e {M e )CVD meet transversally. Moreover, let Z* e and P*. c Al e be the ( m+n )- and m-dimensional 


manifolds obtained by flowing initial conditions on Z and P e forward in (2.4 1 . Finally, we denote 

b y 

J/ £ := y W&tp) c W"(AU 

the union of unstable fibers in ( |2.4[ ) with base points in P*. Note that J/ £ is locally invariant by 
Fenichel’s theory lf2l] Theorem 9.1]. By combining|(E2)|with the Exchange Lemma lf49l Theorem 


49 























2.3], there exists an (m + «)-dimensional submanifold Zi, £ of Z* E and an e-independent neighbor¬ 
hood A c D of ( uq , -fJ(uo), 0,0) such that the Hausdorff distance between An J/ £ and Z i, e is 
0 (e~ c "' f: ). Moreover, trajectories crossing Z\, £ remain in D during the excursion from Z to Zi, £ - 


The singular limit of J/ £ 

First, recall that VI. intersects 7~_ transversely at (uq,- ff(uo))- Second, the unstable manifold 


Wit (Mo) in (2.51 intersects ker(7 - R ) transversely in an m -dimensional manifold <Sq := ]0h(O, u) 


€ U h ] by assumption (El) The a-limit set of So equals the touch-down manifold 7T in At 


by Remark 2.9 We now put these two items together and conclude that the (m + n)-dimensional 


union, 


d/o := U W o(T) c Wq(M,>), 


of unstable fibers in (2.51 with base points in V* () intersects the (m + n )-dimensional manifold 
ker(7 - R) transversally in the point 0h(O, «o). 

Obtaining the periodic orbit 

By Fenichel theory Ell Theorem 9.1] the manifolds J/ £ and J/o have Hausdorff distance (3(e) in 
a neighborhood of the intersection point 0h(O, uf). Therefore, provided s > 0 is sufficiently small, 
J/ £ intersects ker(7-7?) transversally in some point 7\ e , which lies (3(e)-close to 0h(O, mq)- Denote 


by 0h, e (x) the solution to (2.4» with initial condition 0h,e(O) = 7\ £ 


Since 0h(O, «o) is homoclinic to (uo,-J~(uo),0, 0) e At, there exists xq < 0 such that <fih(xo, uo) 
is contained in the neighborhood A of (i<o, -fJ(uo), 0,0). Hence, since 0h,c(O) is (3(e)-close to 
<Ah(0, uo) and xo is independent of s, one derives via Gronwall type estimates that 0h,e(*o) is con¬ 
tained in A n J/ £ . Recall that the outcome of the Exchange Lemma is that J/ £ has Hausdorff 
distance (9(e~ c °A) from Zi, £ in the neighborhood A of <ph, £ ( x o)- 


Denote by Z\ e the (m + «)-dimensional manifold obtained by flowing Z\. E forward in (2.4). 
Since xo is e-independent, we infer, again via Gronwall type estimates, that the Hausdorff dis- 
tance between J/ £ and Z\ £ is 0(e~ Co ^ £ ) in a neighborhood of fh.JO). Therefore, Z\ intersects 
ker(7 - R) transversally in some point P ps , which is (9(e)-close to 0h(O, uf). The solution <p v A x ) 
with initial condition 0 P;£ (O) = P p , £ is the desired periodic orbit. Indeed, 0 Pi£ (x) crosses ker(7 - 7?) 
at x = 0 and at some point x = -L £ < 0, since A p e is contained in Z* E - In particular, we have the 
desired symmetry properties: 0 p , £ (x) = R(t> p _ F (-x) and f VE (L E - x) = Rc/> pE (L E + x) for x e R. 


Checking conditions 1-3 of (S3) 


First, assertion (S3) -2 is immediate by the fact that 0h(O, no) lies (9(e)-close to <'/> p , e (0). Assertion 


(S3) -3 requires more work. By a Gronwall type estimate, we establish 


\\<ph(x. «o) - <PpA x )\\ = 


(A. 2) 


for x e [xo,0]. By the construction in the latter paragraph, we may without loss of generality 
assume that 0 p , £ (xo) is contained in Zi, £ - So, r/> p . £ (x) is in 1 ) for x 6 [-L £ , xo]. Denote by 
a p £ (x), bp.Jx) and c Pi£ (x) the components of the solution r/> p , £ (x) in Fenichel coordinates (A. 1 1 . By 
OT1 Corollary 1] there exists e-independent constants C a , C b ,po > 0 such that 

IIM*)II ^ C b e~^ L C ||a p , e (x)|| < C u e^ (x - x °\ (A.3) 

for x e [-L £ ,xq]. Recall that At £ has Hausdorff distance 0(e~ c °^ £ ) to Ato c At. Hence, using 


estimates (A.2i and (A.31 and the symmetry of </» p£ , estimate (2.91 in assertion (S3)-3 follows. 
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By (A.3l the norm of the tensor a p Jx) 0 b pE (x) is 0(e~ fJnL,: ) on [~L s ,x o]. Therefore, using 
Gronwall type estimates once again, there exists a solution (0,0, c S;£ (x)) in the invariant manifold 
Al £ c {a = b = 0} satisfying d x c = eK(c,e), which is (9(e -W)L£ )-close to c PiE (x) for x e [-L e ,x o]. 
Moreover, there exists a solution c s ,o(x) to d x c = K{c, 0), which satisfies ||c p e(x) - c Si o(e~ 1 x)|| - 
(9(e) for x e [-L E ,x q]. Recall that system d x c = K(c, 0) corresponds to the slow reduced system 
\2.6\ . Thus, we have established the last estimate ( 2. 10[ ) in |(S3)| -3 and assertion (S3) T. Note that 
the solution c S) o(n) corresponds to the solution if/ s (x) to (2.6) in (E2) □ 


B Prerequisites for the proofs of the main results 

In this section we treat some prerequisites needed for the proofs of our main results. We will 
provide the proofs of statements that we could not find in literature. 

B.l Bounded growth estimate 

We start with the notion of bounded growth, which is due to Coppel 191 . 

Definition B.l. Let n e Z>o, /c Ran interval and A e C(J, Mat nx „(C)). Denote by T(x,y ) the 
evolution operator of system. 


d x tp = A(x)ip, <peC n , 


(B.l) 


Equation (B.l I has bounded growth on J with constants p, K > 0, if for all x ,y e 7 it holds 

||r(jc,y)|| < K^ x ~ y \ 


Clearly, if A(x) is bounded on R., system (B. 1 1 has bounded growth on R with constants p = ||A|| 
and K = 1 due to Gronwall’s inequality. In our spectral analysis, we often want to compare a linear 
system with its perturbation. This requires the following consequence of Gronwall’s inequality for 
linear systems. 


Lemma B.2. Let n e Z>o, a.b e R, a < h and A,B e C([n, b], Mat, lx «(C)). Denote by T(x,y ) the 

(B.2) 


evolution operator of system (B.l). Similarly, denote by T'(x,y ) the evolution of system, 

d x (p = B(x)ip, ip e C”. 


Suppose (|/i. 1\ has bounded growth on \a, ft] with constants p, K > 0. It holds 


II T(x,y) - T'(x,y )|| < K J ||A(z) - 5(z)||*exp|q(ft _ a ) + K J 
for all x,y € [a, ft]. 

Proof. This follows from the proof of ll4Tl Lemma 1], 


||A(z) - B(z)\\dz\, 


B.2 Exponential dichotomies 

Exponential dichotomies are an important tool in studying spectral properties of differential equa¬ 
tions. They enable us to track solutions in the linear stability problem by separating the solution 
space in solutions that either decay exponentially in forward time or else in backward time. 
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Definition B.3. Let n e 


/c Ran interval and A e C{J, Mat„ x „(C)). Denote by T(x,y ) the 
evolution operator of (B.l). Equation (B.l i has an exponential dichotomy on J with constants 
K,p > 0 and projections P(x) : C" —» C", x e J if for all x,y e J it holds 


• P(x)T(x,y) = nx,y)P(y)-, 

• ||r(jc,y)/ 5 (y)|| < Ke~^ {x ~ y ^ for x > y\ 

• ||r(x,y)(/ - P(y))\\ < Ke-^ y - x) for y > x. 


Let P(x), x € / be the family of projections corresponding to an exponential dichotomy on J. For 
each x € 7, the range P(x)[C"] and the kernel ker(P(x)) are called the stable and unstable sub¬ 
spaces at x, respectively. These spaces are often denoted by E s (x) and E“(x), leaving the projection 
P(x) implicit. Similarly, we abbreviate T s (x,y) = T(x,y)P(y) and T u (x,y ) = T(x,y){I - P(y)). 

Below we give a short overview of the properties of exponential dichotomies that we need for 
our spectral analysis. For an extensive introduction on dichotomies the reader is referred to the 
work of Coppel 0. A generalization of the concept of exponential dichotomies is the notion of 
exponential separation, which is treated in a paper of Palmer |f39l . 


B.2.1 Sufficient criteria 

Clearly, an autonomous linear system with hyperbolic coefficient matrix has an exponential di¬ 
chotomy on R. This result can be extended to non-autonomous linear systems under the hypothesis 
that the coefficient matrix changes sufficiently slow and remains bounded. 

Proposition B.4. Let n e Z>o, J c Ran interval and A 6 C ] (J, Mat„ X /j(C)) such that 

1. There exists a > 0 such that for each x £ J the matrix A(x) has (coun ted with algebraic 
multiplicity) k eigenvalues with real part < —a and n — k eigenvalues with real part > a. 

2. There exists M > 0, which bounds A on J. 

There exists 6 > 0, depending only on a and M, such that, if for all x € J we have ||A'(x)|| < 5, then 
( B.l\ has an exponential dichotomy with constants K,p > 0 and projections P(x), x € J. Here, 
we have p = \a and K depends only on M and a. Moreover, there exists a fundamental matrix 
solution X(x) oft \ll I\ such that 

P(x) = X(x) | ^ o) X " 1W ’ XeJ - 

Proof. This is the content of 0 Proposition 6.1]. □ 


B.2.2 Extending and pasting 

Once one puts a linear stability problem in the framework of exponential dichotomies, a great 
technical toolbox becomes available. Namely, on the one hand, there is a scala of constructions 
available to extend the interval of the dichotomy. On the other hand, exponential dichotomies 
persist under small perturbations of the equation. Therefore, these techniques enable us to de¬ 
compose our linear stability problem into simpler subproblems. In this section, we treat extension 
and ‘pasting’ of exponential dichotomies. In the next section |B.2.3 we consider the persistence of 
exponential dichotomies under small perturbations. 
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Lemma B.5 (Extension Lemma). Let n £ Z>o, Ji c J\ c R internals and A e C(.J \, Mat raxn (C)). 


Suppose equation (B.l) has an exponential dichotomy on .L with constants Ki,p > 0 and projec¬ 


tions P 2 (.x), x £ J i- Then, if the length of J \ \ J 2 is finite, system (B.l) has an exponential dichotomy 
on J\ with constants K\,p > 0 and projections P\ (x), x £ J\. The constant K\ satisfies 


7^= 7^2 exp [ j (2\\A(x)\\ + p)dx 


Moreover, we have P\{x) = Pfipdtfor all x £ 72- 


Proof. This is proven in J9[ p. 13]. □ 

The next result shows that, if a periodic equation admits an exponential dichotomy on a sufficiently 
large interval, then it has an exponential dichotomy on the whole line. 


Lemma B.6 (Periodic Extension Lemma). Let n £ Z>o, T > 0 and A e C(R, Mat, lx; ;(C)). Suppose 
that A is T-periodic and that equation (B.l I has an exponential dichotomy on an interval J of 
length 2T with constants K,p > 0. Let M > sup ce? ||A(x)|| and h = /r _1 (sinh _1 (4) + log(TQ). 


If T > 0 is so large that T > 2h, then equation < \B. /[ ) has an exponential dichotomy on K. with 
constants K\,p\ > 0. We have p\ = h~ l log 3 and K\ depends only on M, K andp. 


Proof. This is the content of fid] Theorem 1]. □ 

Exponential dichotomies on two connected intervals can be pasted together as long as their stable 
and unstable spaces at the connection point are complementary. 


Lemma B.7 (Pasting Lemma). Let n € Z>o, 7 c R an inten’al and A 6 C(7,Mat„ x „(C)). Let 
Ji, J 2 be two intervals such that their union equals J and max J\ - b = min J 2 for some be R. 
Suppose equation (B.l) has exponential dichotomies on both 7i and J 2 with constants K,p > 0 
and projections P\(x), x £ J\ and P 2 (x), x £ J 2 , respectively. 


If E'fib) = kci( P\(b)) and E^(b) = 7S (/?) IC” | are complementary, then (B.l) has an exponen¬ 
tial dichotomy on J with constants K\,p > 0. Here, K\ depends only on K and ||P||, where P is the 
projection on Efib) along Efib). 


Proof. Let X(.r) be the fundamental matrix of (|B. 1|) satisfying X(b) = 7. Define P(x) = X(x)PX 1 (x) 


for x £ 7, where P is the projection on Efib) along E"(b). Observe that P = P(b) has the same 
range as P 2 {b) and the same kernel as P\(b). Now, the exposition in (9] pp. 16-17] shows that 


(B.l 1 has exponential dichotomies on J\ and on J 2 with constants K\,p > 0 and projections P(x) 
for x £ J\ and x £ J 2 , respectively. We have K\ = K + K 2 \\P\\ + K 3 . To conclude the proof we need 
to show that the dichotomy estimates remain true on the union 7 = 7i U 72- Indeed, take x e J 2 
and y £ 7]. We estimate 


||7Xx,y)P(y)|| < \\T(x,b)P 2 (b)\\\\P\\\\P 1 (b)T(b,y)\\ < ^ 2 ||P||«-^>, 


where we have used P 2 (b)P = P and PP\(b) = P. Similarly, one estimates \\T(y,x)(I - P(x))|| < 
A' 2 ||/j|e~ / ' ( * _v ) for x £ 72 and y £ 7|. □ 
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B.2.3 Roughness 


Exponential dichotomies are in particular useful to study the spectral properties of perturbed dif¬ 
ferential equations, since they persist under small perturbations of the equation. This property is 
often referred to as roughness or robustness. We start with a general roughness result for exponen¬ 
tial dichotomies on arbitrary intervals. 


Proposition B.8 (Roughness on arbitrary intervals). Let n e Z>o ,6 > 0, / c I an internal and 
A,B e C(J, Mat„ x „(C)) such that 


sup \\A(x) - R(x)|| < S. 

xej 


Suppose (B.l) has an exponential dichotomy on J with constants K,p> 0 and projections P(x) : = 
X(x)PX (xj, x £ J, where X(x) is a fundamental matrix of \B. /[ ) and P an orthogonal projection. 
If we have 


6 < 


36K 5 ’ 


then equation ( |ff.2[ ) has an exponential dichotomy on J with constants K\,p\ >0 and projections 
P](x), x e J satisfying p\ - p- 6 K 3 6, K\ - 12 K 3 and for all x £ J 


\\P(x)-Pi(x)\\< 


144 K 6 6 
P 


Proof. This is the content of l9l Proposition 5.1]. □ 

The latter persistence result can be simplified significantly in the case / = R. 


Proposition B.9 (Roughness on R). Let n e Z>o,d > 0 and A,B € C(R, Mat„ xn (C)) such that 

sup ||A(jc) - B(x)|| < 6. 

X€R 


Suppose equation (B.l I has an exponential dichotomy on R with constants K,p > 0. If we have 

P 


5< 


4 K 2 ' 


then equation (B.2 1 has an exponential dichotomy on R with constants K\,p\ > 0. Here, we have 
p\ = p — 2 K6 and K\ depends on K only. 


Proof. This can be found in (9l pp. 34-35]. □ 

If an equation has an exponential dichotomy on R, then it admits no non-trivial bounded solutions. 
It is possible to achieve persistence of the latter fact under milder conditions than those stated in 
Proposition |B.9| 


Proposition B.10 (Persistence of no non-trivial bounded solutions). Let n e Z>o and A,B € 



we have 


\\Ti(x,y)-T 2 (x,y)\\<l, 


then (B.2) admits no non-trivial bounded solutions. 


Proof. This is the content of Il42l Theorem 1]. 


□ 
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B.2.4 Uniqueness 


We emphasize that exponential dichotomies on an interval / c R are in general not unique. For 
instance, if J = R>o, then the stable subspace E s ( 0) is uniquely determined, whereas the unstable 
subspace E u ( 0) can be chosen to be any complement of E s (0). Flowever, given two exponential 
dichotomies it is possible to estimate the ‘gap’ between the stable subspaces and unstable sub¬ 
spaces. 


Lemma B.ll. Let n e Z>o, a,b e R with a < b and A,B e C([a, b\, Mat„ xn (C)). Suppose 
equations (B.l I and (B.2 1 have exponential dichotomies on [ a , b ] with constants K 1 , 2 ,/q ,2 > 0 and 
projections P\Ax), x £ Denote by T \ 2 (x, y) the evolution operators of systems (B.l) and (B.2). 
Let 6 > 0 such that 


\\T 1 (a,b)-T 2 (a,b)\\<6. 

Then, for every v e Efa) = P\(a)\C" |, there exists w e Efa) = P 2 {a)\C n ] such that 

||v - w\\ < (6 + K 2 e-^ b - a) )K ie - Mb - a) \\v\\. 


(B.3) 


Similarly, for every v £ E“fb) = kcv(P \(b)), there exists w £ E“(b) = kcr( P 2 (b)) such that (B.3) 
holds true. 


Proof. Let v e Efa) and consider w = T 2 (a, b)P 2 (b)T \ (b, a)v e Efa). We estimate 

\\T 2 (a,b)P 2 (b)Tfb,a)v - v|| < [|| T 2 (a,b) - Tfa,b)\\ + \\T 2 (a,b)(I - P 2 (b ))||] ||ri(&,a)v|| 

< (6 + K 2 e-^ b - a) )K ie -^ b - a) \\v\\. 

The other statement is proven in an analogous way. □ 


B.2.5 Inhomogeneous problems 

When comparing solutions to inhomogeneous problems using Gronwall’s inequality, one often ob¬ 
tains sharp estimates on finite intervals only, since the solutions to the corresponding homogeneous 
problems will grow in general exponentially. Exponential dichotomies prove to be important tools 
to compare bounded solutions to inhomogeneous problems on the whole real line. This is the 
content of the following technical result. 


Proposition B.12. Let n € Z>o,/, g £ C(R, C") bounded and A,B e C(R,Mat„ x „(C)). Suppose 
equation \B.l ) has an exponential dichotomy on R with constants K,p > 0. Then the inhomoge¬ 
neous problem, 


dxT = M*)<P + fix), tfi £ 


(B.4) 


has a unique bounded solution (fix). Furthermore, suppose in addition that A and B are bounded. 
Let a,b £ R with a < b. Then, for any bounded solution ijj(x) to the inhomogeneous problem, 

8 x cj = B(x)a> + g(x), cj £ C”, (B.5) 

we estimate for x £ [ a , b] 

Mx) - <K*)II < - (e^ {x ~ a> + e-* b ~ x) ) (M||A - fi|| + ||/ - *||) 

p v ' 

2k ( ) (B-6) 

+ — 11/11 sup ||A(z) - fi(z)|| + sup ||/(z) - *(z)|| . 

P V Z€[a,b] ZE[a,b] ) 
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Proof. Denote by T{x,y) the evolution of system ( B.l[ ). By (9J Proposition 8.2] system (B.4i has 
a unique bounded solution given by 


X X poo 

T\x, z)f(z)dz - I T u (x, z)f(z)dz, 

CO kJ x 


x e. 


Now, let A and B be bounded and ip a bounded solution to (B.51. Note that w: 
w(x) = tfi(x) - f{x) is a bounded solution to the inhomogeneous equation, 

d x w = A(x)w + h(x). 


l' 1 defined by 


where the inhomogeneity h: R —» C" given by h(x) = (A(x) - B(x))fix) + fix) - g{x) is bounded 
on R by hypothesis. By applying |9] Proposition 8.2] once again we deduce that w{x) is given by 


-I 


w(x ) = T\x, z)h(z)dz 


s>oo 

T u (x, 

Jx 


z)h{z)dz, x 6 


(B.7) 


Now, let a, b € R with a < b. Estimate (B. 61 for v e \a, b \ is achieved by splitting both integrals 
in expression S3 into two parts. The first integral is split in integrals over (-co, a) and over 
(a, x). Similarly, the second integral is split in integrals over {x, b) and over (b, oo). This yields 
four integrals, which can be estimated separately in order to obtain estimate (B. 61 . □ 


B.3 The minimal opening between subspaces 

The minimal opening |[24l Section 13.3] is a quantity measuring the ‘gap’ between two subspaces. 


Definition B.13. Let n e Z>o- The minimal opening between two non-trivial subspaces At and N 
of C" is given by 


p(M, AO = inf{11 -T - y|| : x e M,y £ N, max(||.r||, ||y||) = 1}. 


The minimal opening has the useful property that the norm of the projection on At along U can 
be bounded in terms of rj( At, N). This norm estimate is essential for the application of the pasting 
Lemma B.7 in our spectral analysis. 


Proposition B.14. Let n € Z>o. The following assertions hold true. 

1. If P is a non-trivial projection on C n , then it holds 

||P|| < ---. 

77(P[C n ],ker(P)) 

2. For non-trivial subspaces At and N of C" it holds ?/(At, N) T 0 if and only ifM D N = (0). 

3. Let At 12 and N\ i be non-trivial subspaces ofC n . Suppose that there exists 0 < 6 < 1 such 
that for each v e At,- there exists aw € At, such that ||v - w\\ < d||v||/or i = 1,2. Then, we 
have the estimate 


r](N \, No) < r)(M ] , M 2 ) + 46. 

4. Let LI c C be open and connected. Suppose Mi A) and NiA) are continuous families of 
subspaces on Q, i.e. there exist continuous families ofprojections PPu ■ Mat„ x „(C) 
such that P;vi(/l)[C' ! ] = Mi A) and Pyv(/t)[C!' ! ] = N{A) for A e LI. Then, the map A i—> 
rj (yVt(/l), A ({A)) is also continuous on LI. 
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Proof. The first two assertions are derived in (2J[ p. 396] and mi Proposition 13.2.1], re¬ 
spectively. For the third assertion take e > 0. There exists vi e M\ and v 2 £ M 2 with 
max(||vi||, ||v 2 ||) = 1 such that ||vi - vi\\ < t]^M\,M 2 ) + c. Without loss of generality we may 
assume ||vi|| = 1. By hypothesis there exists w 1 £ A/j such that ||vi - will < 8. Because we have 
8 < 1, we can normalize w\ and define z\ := pp. One readily estimates ||vi — ziII < 28. Similarly, 
there exists w 2 £ N 2 such that ||v 2 - W 2 II < 8. In the case ||w 2 l| > 1, take Z 2 pp. One easily 
verifies ||v 2 - Z 2 II < 28. In the case HwtII < 1, we just take Z 2 ■- w 2 - Finally, we estimate 

t](Ni,N 2 ) < ||zi -Z 2 II < llvi - v 2 || + ||vi — zill + ||v 2 — Z 2 II < ti(M\,M 2 ) + 48 + e. 

Since e is arbitrarily chosen, the second assertion follows. Finally, for the fourth assertion let 
PA) and P,\AA) be continuous families of projections on Q with ranges yYl(d) and N(A), respec¬ 
tively. With the aid of identities (13.1.4), (13.2.5) and (13.2.7) in If24l we derive for Ao £ Q 

1 7](M(A),JV(A)) - 7j(M(Ao),JV(Ao))l < V2(|| \P M (A) - P M (^o)ll + II^U) - ^Uo)ll) ■ 

This shows that A —> r](M(A), N(A)) is continuous on fl □ 


C Proofs of technicalities 


Lemma C.l. The map Q, defined in (3.8), is bounded at 00 . 


Proof. The coordinate change (v, q) (v, Vdw) puts system (3.5) into the form, 


D 2 d x v = s/~Aw 

q w _ I d v G(uo,Vh(x),0) 

V2 


+ 


vr)v ■ ”• 


W £ 


(C.l) 


If A > 0 is sufficiently large, system (C.l l has by Proposition B.9 an exponential dichotomy 
on K. with constants K\,p(A) > 0, where K\ is independent of A and p(A) = ^|| VO 2 II 1 Vd. 


Therefore, system (3.5) has an exponential dichotomy on R with constants K 2 (A),p(A) > 0, where 
K 2 = Vdk'j and K\ is independent of A. Note that is d-independent. Moreover, since G 
vanishes at v = 0 by |(S1)| ^ 21,0 is exponentially localized by Remark 2.4 Combining these fact 
with Proposition |B . 12 1 yields that, for d > 0 sufficiently large, there exists K 2 ,p 3 > 0, independent 
of d, such that \\Xi n (x, d)|| < K^e~ M3 ^. This concludes the proof. □ 

Lemma C.2. Let m = 1. The trace t(A) = Tr(T(d)7~ s (2Lo, 0, d)) diverges to 00 as .1 ^ 00 . 


■ P e 


Proof. Take d > 0. Consider system, 

dpi = sfXp 

dx p = ^-(Ms(^),0,0) + Vdjn 

with evolution T s \(x,y, d). Denote by T s2 {x,y, d) the evolution of the autonomous system, 

I d x u = sfXp 


(C.2) 


\ dsp = Vd(i 


u, p e 


(C.3) 


Let M > 0 be a bound of 0) on [0,2Lo], Syste m (|C.3| ) has bounded growth with 

constants K = 1 and p = Vd. Therefore, it holds by Proposition B.2 for d > 0 sufficiently large 


l|7~ s i (2L 0 ,0, d) - T s2 (2Lq, 0, d)|| < 2ML « gM g 2VdL 0 

Vd 


(C-4) 
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On the other hand, system ( |3.7| ) is equivalent to system ( |C.2[ ) upon performing a coordinate change. 
Indeed, take C,i '■= (q ^). It holds 


C A T S] (2L 0 , 0, A)C~/ = T S (2L 0 , 0, A), (C.5) 

We approximate t(A) for A > 0 sufficiently large with the aid of Lemma |CTj (|C.4[) and (C.5|) 


t(A) = tr 


1 0 

0A) , 


I Ts\(2L(), 0, A) 


T s2 (2Lo,0,A)\ + 0\^e 2 ' rAh 


= e 


1 ° 

gW) 1 

VI ) / 

2 VIio + e 2VlLo + _J_^(q)(g2VlLo _ g-2VlLo) + Q (<? 2 ^ j 


_ 2 Vlfo 


= e 


1 +o 


VI 


In the latter approximation we have used explicit expressions of the evolution T s2 {x,y,A) of the 
linear autonomous system (|C.3[). We conclude t{A) —> oo as A — > oo. □ 


Proof of Proposition [67TT} 2. Let A > 0. Putting y = \[Ax and p = x] AD \ r rescales system (3.7 1 
into 


u = r 


y/Didyr = 


1/2 y),0,0) 


+ l\u ’ 


u, r e 


Denote by 7~ v i (>’, z. A) the evolution of system (C. 61 . It holds 

C A rAX)T s A2^AL {h (),X)C- x ' = T(/t)7 _ s (2L 0 ,0, A), 


(C.6) 


(C.7) 


with C x := (o vJjt) and Ti(d) := ( | vizvl'gw) / )• System (C.6) is close to, 


\[U\dyU = 

VDTd'i 


r = u 


u, r e 


(C.8) 


Clearly, (C .81 has an exponential dichotomy on R. with constants K\ — I and jo = ||Z>i || l/2 . The 
corresponding rank m projections P\ = I ( / ) arc independent of x, since (C .81 is autonomous. 

Let M \ > 0 be a bound of d n Hi(u s (-), 0,0) on [0, 2Lq |. So, by roughness (Proposition B .81 system 
(C .61 has, for A > 0 sufficiently large, an exponential dichotomy on [0,2 VdLo] with constants 
K 2 ,p 2 > 0, independent of A. The corresponding projections P 2 (x, A), x e [0,2 VaLq] satisfy 


\\P 2 (x,A)-Pi\\ < 


144 M 


Pi A 


(C-9) 


Now, choose bases B { ' 6 Mato, 


of Pi[C 2m ] = P;|C" ! ] andker(Pi) = S“[C m ]. Define 


B s 2 (A) = P 2 ( 0, A)B\ and B“(A) = (I - P 2 (2 VdL 0 , A))B\. By estimate (jC9| it holds 

T 


P^(T) -B u { s \\ = 0 1 

Consider the invertible matrix, 

'H(A) ;= (t~ s i (0,2 VaL 0 , A)B u 2 (A), B s 2 (A )). 


(C.10) 


By Lemma C.l and (C.10) we have 

(T 1 (T)r; s i(2 VaLo, 0,4) - y) 'H(A) = (B\, -yB\) + O |-J= 
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Taking determinants in the previous expression gives by (C.7 1 


6,. 0 (A, y) detc H{A)) = det (T 1 (A)7“, 1 (2 VdL 0 , 0 , 4) - y) detCK(d)) = (-y) m det (b\, + <9 

By construction det(B", fij) is non-zero and independent of A. Combining this with dctbH(A)) + 0 
yields that £. v ,o(9, y) must be non-trivial for A > 0 sufficiently large. □ 
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